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ABSTRACT 

The performance of the Learning Fuzzy Logic Control System (LFLCS), devel- 
oped in this thesis, has been evaluated. The Learning Fuzzy Logic Controller (LFLC) 
learns to control the motor by learning the set of teaching values that are generated by 
a classical PI controller. It is assumed that the classical PI controller is tuned to mini- 
mize the error of a position control system of the D.C. motor. The Learning Fuzzy 
Logic Controller developed in this thesis is a multi-input single-output network. 
Training of the Learning Fuzzy Logic Controller is implemented off-line. Upon com- 
pletion of the training process (using Supervised Learning, and Unsupervised Learn- 
ing), the LFLC replaces the classical PI controller. In this thesis, a closed loop 
position control system of a D.C. motor using the LFLC is implemented. The primary 
focus is on the learning capabilities of the Learning Fuzzy Logic Controller. The 
learning includes symbolic representation of the Input Linguistic Nodes set and Out- 
put Linguistic Nodes set In addition, we investigate the knowledge-based representa- 
tion for the network. As part of the design process, we implement a digital computer 
simulation of the LFLCS. The computer simulation program is written in “C” com- 
puter language, and it is implemented in DOS platform. The LFLCS, designed in this 



thesis, has been developed on a IBM compatible 486-DX2 66 computer. First, the per- 
formance of the Learning Fuzzy Logic Controller is evaluated by comparing the angu- 
lar shaft position of the D. C motor controlled by a conventional PI controller and that 
controlled by the LFLC. Second, the symbolic representation of the LFLC and the 
knowledge-based representation for the network are investigated by observing the 
parameters of the Fuzzy Logic membership functions and the links at each layer of the 
LFLC. While there are some limitations of application with this approach, the result 
of the simulation shows that the LFLC is able to control the angular shaft position of 
the D.C.motor. Furthermore, the LFLC has better performance in rise time, settling 
time and steady state error than to the conventional PI controller. 

This abstract accurately represents the content of the candidate’s thesis. I recommend 
its publication. 
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Chapter 1 


1.1 Introduction 

In the design of Analog and Digital Control systems, a dynamic representation of 
the system is required. In design engineering, this information is usually not known as 
a priori. Furthermore, in order to minimize the errors between the output and the 
input, modem control theory requires that an accurate model of the system be avail- 
able. These requirements limit the application of modem control theory in many areas. 
It is important that controllers be developed that do not have such stringent require- 
ments. The goal of this research is then to explore an alternative controller and to eval- 
uate its performance. In this thesis, the concepts of Fuzzy Logic Control and Neural 
Network Learning are combined to design a controller for an unknown plant. The 
developed system is referred to as a Learning Fuzzy Logic Control System (LFLCS). 

Whereas Unsupervised Learning Neural Network is used to set up the initial struc- 
ture for the network controller, Supervised Learning Neural Network is used to adjust 
parameters of the network controller to minimize the output error. This research 
addresses the learning capability of the Learning Fuzzy Logic Controller (LFLC), and 
its knowledge representation in symbolic terms. A complete system diagram is shown 
in Fig.l. The LFLCS requires that the input signals and the teaching signal be avail- 
able for training, and that a good logical structure be set up before the training takes 
place. The LFLC is rather different from a conventional controller; this difference is 
explained and illustrated. 
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1.2 Previous Research 

Modem control theory has proved to be very useful in areas where systems are 
well defined either deterministically or stochastically. However, many control systems 
involve human-judgement interaction as part of the control system. Human involve- 
ment often provides an adequate controller because the mind of an operator usually 
provides a model of the process which is just accurate enough to carry out the task at 
hand. On the contrary an automatic controller has no way of observing the essential 
features of a particular process. A human is usually capable of learning through expe- 
rience, which decreases the need for a precise model of the system. Thus, modelling 
the human decision making process is essential in control system design. The knowl- 
edge of the control process in the human mind is captured in the fuzzy system design 
approach. 

In 1965, Zadeh [31] introduced and developed the concept of fuzzy set theory. 
Since its introduction, fuzzy logic has been successfully applied in many control sys- 
tem applications; for example, see references [1], [29], [6], [17], [16], [20] for some 
good illustrations. An excellent overview of fuzzy logic applications in control engi- 
neering is given by Langari and Berenji [15]. In reference [4], Chuen-Tsai Sun imple- 
mented the fuzzy IF-THEN rule base to identify the structure and the parameters of a 
network such that a desired input-output mapping is achieved. However, a fuzzy logic 
controller requires that the control strategy be obtained as a fuzzy set term. This limits 
somewhat the usefulness of the concept of a fuzzy controller. 

The Self-Organizing (unsupervised learning) approach was presented by Zhang 
and Edmunds [32], and Linkens and Hasnain [20]. They proposed a fuzzy logic con- 
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trailer that is able to cluster (self-organize) the input data without any prior knowledge 
of the data, and the network automatically sets up the parameters for each membership 
of the network. Furthermore, this network has a learning algorithm and is capable of 
generating and modifying control rules based on an evaluation of the system perfor- 
mance. The generation and modification of the control rules is achieved by assigning 
a credit or reward value (weight) to the individual control action that makes a major 
contribution to the current performance. This is an excellent control strategy for a sys- 
tem when the operator control strategy is not availble. 

The concept of machine learning was introduced many years ago in an effort to 
achieve human like performance in the fields of speech and image recognition for 
handicapped people. An extensive research effort to simulate the process of intelligent 
human learning using an Artificial Neural Net using the neural network is a major 
effort in the field of computer sciences [2], [9], [11], [12], [24], [30], [26]. Also, Self- 
Organized (unsupervised learning) is one of the learning methods which is used in 
speech recognition. Supervised learning is used in many fields where teaching data 
are obtainable from an applicable sources [21], [22], [10], [27], [8], [25]. 

Regardless of the name these models attempt to achieve a good performance via a 
dense interconnection of simple computational elements. In this respect, an artificial 
neural net structure is based on present measurements of biological systems. Neural 
net models have the greatest potential in application such as speech, image recogni- 
tion, and control systems where many of the hypotheses require a parallel, high com- 
putational rate, and the better systems result from observing inadequate human 
performance. It is hoped that the potential benefits of neural nets will extend beyond 
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the high computation rates provided by present massive parallel system. Neural nets 
typically provide a greater degree of robustness or fault tolerance than the von Neu- 
mann sequential computers of today because there are many more processing nodes 
possible, each with primarily local connections. The presence of a few imperfect 
nodes or links thus need not impair the overall performance significantly. Most neural 
net algorithms also adapt connection weights in order to improve the performance 
based on current results. Work on artificial neural net models has a long history. 
Development of detailed mathematical models began more than 45 years ago with the 
work of McCulloch and Pitts [23]. Lin and Lee presented in [18] a two-phase learning 
fuzzy logic network which consisted of both unsupervised learning and supervised 
learning, and in [19] they developed a reinforcement neural network-based fuzzy logic 
control system. 

In this thesis, following the approach of Lin and Lee [18], the Learning Fuzzy 
Logic Control System (LFLCS) is proposed in which the learning capabilities of neu- 
ral networks are utilized. The system learns by adjusting the parameters of the neural 
network using training data. The learning schemes of a Learning Fuzzy Logic Con- 
troller (LFLC) combine both the unsupervised (self-organized) and the supervised 
learning. 

The neural network structure, implementing the LFLC, is given in Fig.2. This net- 
work has one output and two inputs. The input signals x\, and x 2 to the LFLC are the 
feedback signal and the error signal of the control system, respectively, as shown in 
Fig.l. The teaching pattern used by the LFLC is the control signal generated by a con- 
ventional controller. 
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Fig.l. Learning Fuzzy Logic Control System 



1.3 Structure of the LFLC: An Overview 

A Neural Network (or connectionist network) is a highly parallel connected net- 
work that attempts to model the learning ability of the human brain. The intelligence 
of the network is represented by the weights that connect nodes at one layer to the 
nodes of the next layer of the network. Fuzzy Logic Control is a knowledge-based 
control strategy. This strategy can be employed given a sufficiently accurate control 
law which is not unreasonably complex. The Neural Network learns by the tuning of 
system parameters using training data. The learning schemes of the LFLC combine 
both the unsupervised (self-organized) learning and the supervised learning. A layout 
of a simple network is found in Fig. 2. This particular network has two inputs and one 
output, i.e., it belongs to a Multiple Input Single Output (MISO) class of networks. As 
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shown in Fig.l, the input signal (xj) to the LFLC is the feedback signal of the control 
system, and the input signal (x 2 ) to the LFLC is the error signal of the control system. 
The LFLC learns to recognize a set of data. This set of data is called the teaching pat- 
tern which is generated by the digital Proportional-Integral (PI) controller in Fig. 1. 

The Neural Network is used as a learning mechanism and the Fuzzy Logic Control 
algorithm is actually controlling the plant The LFLC is a two phase learning network. 
The first phase is an unsupervised learning phase, and the second is supervised learn- 
ing phase. The LFLC has a total of five layers. The first layer is the Input Linguistic 
Nodes layer which in this particular implementation contains two nodes. These nodes 
represent the input data sets as symbolic terms. A set of five nodes is set up in the sec- 
ond layer for each Input Linguistic Node. This is the Input Term Nodes layer. The 
purpose of the Input Term Nodes is to categorize the input data into linguistic terms. 
The third layer is the Rule Nodes layer. Each rule node in the third layer represent s a 
rule of controlling the plant. The number of rule nodes in this layer is equal to the 
product of the numbers of nodes in each set of the Input Term Nodes. Therefore, there 
are 25 rule nodes in our imnlementation. The fourth layer is the Output Term Node 
layer that contains seven nodes. The purpose of the Output Term Node is to categorize 
in linguistic terms the consequences of the fired rules. The fifth layer is the Output 
Linguistic Nodes layer that in the case of the application considered contains one out- 
put node. However, a second node is used to train the Output Linguistic Node in the 
unsupervised learning phase. The Output Linguistic Node constitutes the output of the 
LFLC. The data are randomly presented to the LFLC in the learning phases. How- 
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ever, before the data are presented to the network, they go through a normalization 
process, which is explained in detail in the next section. 


7 



Chapter 2 


2.1 Normalization 

In the normalization process, the data of the error signal (x 2 ) and the feedback sig- 
nal (jcj) are mapped to the range of [-1,1]. The data of the teaching pattern (y t ) are 
mapped to the range of [0,1]. These mappings are accomplished by the data normal- 
ization using the following: 


= 


iW > 

max(\max jc, \j] |, | min X; [/] |) ’ 


i = 1,2 j = l,2,...n 


( 1 ) 


05 y/U} 

* ’ 2 • max ( | max yj ,|min y \ ) ’ 


j = 1,2,. ..n 


( 2 ) 


where n denotes the number of data points. 

Since [jq[/] } contains all data of the linguistic node (xj) that are going to be nor- 
malized by Eq.(l), the {y,[/] } contains all data of the linguistic node (y f ) that will going 
be normalized by Eq(2). The max x t becomes the maximum number of the set {*,[/] } , 
and min x t the minimum number of the set [x,[/]}. In the same manner the max y t is 
the maximum number of the set {)>,[/]}, and min y t the minimum number of the set 
{y,[/]}. The x t is the result of the normalization of the data of the linguistic node (x t ), 
and y t is the result of the normalization of the data of the teaching patterns y t . This 

normalization ensures that the input data and the teaching data are mapped to the 
ranges [-1, 1] and [0,1] respectively. All negative values of the original teaching data 
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are mapped to [0, 0.5), and all positive values of the original teaching data are mapped 
to (0.5, 1). Similarly, all negative values of the original input data are mapped to [-1, 
0), and all positive values of the original input data are mapped to (0, 1]. 


•3 

■) 
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Chapter 3 


3.1 Network Structure 

The LFLC has a total of five layers. These are the following: 

First layer or Input Linguistic Nodes Layer 
Second layer or Input Term Nodes Layer 
Third layer or Rule Nodes Layer 
Fourth layer or Output Term Nodes layer 
Fifth layer or Output Linguistic Nodes Layer. 

Each of these layers is connected by the link between them. The purpose and detail of 
each layer will be explained in the following sections. 

3.2 Input Linguistic Nodes Layer 

The purpose of this layer is to propagate the input data to the Input Term Nodes 
layer; therefore, each input linguistic node is connected to a set of Input Term Nodes in 
the next layer. The output value of the input linguistic node is the same as the input 
value, and is propagated to its own Input Term Nodes set in the next layer. The link 
(w) from first layer to the second layer is unity. 

3.3 Input Term Nodes Layer 

The purpose of the Input Term Node is to represent the value of the Input Linguis- 
tic Node in linguistic terms. For example, the linguistic term for an Input Term Node 
can be Negative Large (NL), Negative Medium (NM), Negative Small (NS), Zero 
(ZN), Positive Small (PS), Positive Medium (PM), or Positive Large (PL). 



Two sets of Input Term Nodes are set up for the LFLC as shown in Fig. 2. Five 
nodes are set up for each Input Term Nodes set; thus, for this particular network, a 
total of 10 nodes are in the Input Term Nodes layer. The input value to each node in 
the Input Term Nodes set number 0 is equal to the product of the output value of the 
Input Linguistic Node number 0 and the link (w) that connects to the Input Term node. 
Each node of the Input Term Nodes set has a membership function. This membership 
function can be based on any activation function (a); however, a Gaussian activation 
function is used in the network considered. This Guassian activation function is 
defined by (3) and (4) below. It is used for all memberships in layer 2, 



v 


where/is the memberships function M J f m... S..Y The parameters m 0 and S i; are the 

V V/ 

center (the mean) and the width (the variance) of the Gaussian function respectively, 
and uf is the input data to the Gaussian function with the superscript 2 used as a refer- 
ence to the second layer. The subscript i is the index of the Input Term Nodes set, and 
the subscript j is the index of the nodes within each Input Term Nodes set. The links 
(w) of this layer are fully connected and they are equal to unity (See Fig, 1.). For our 
particular structure shown in Fig. 2, i = 1, 2 and j = 1, 2, 3, 4, 5. We have u\ = Xj and 

«2 = *2- 
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3.4 Rule Nodes Layer 

This layer contains a set of fuzzy logic rules Rj. For this thesis, a MISO system is 
used. A linguistic variable (x) in a universe of discourse U is characterized by two 

sets: T(x) = {Tj,7? T k } and MU) = {M' x Ml,...M k x } ■ The T(x) is the term 

set of x which is the set of names in linguistic terms of the values of x with each value 
T‘ x being a fuzzy number with membership function M‘ x defined on U. Thus M(x) is a 
semantic rule for associating each value with its meaning. 

For example, if x indicates voltage, then T(x) may be in the set of {NL, NS, ZN, 
PS, PL}. In this thesis, five Input Term Nodes are set up for each Input Linguistic 
Node, and seven Output Term Nodes are set up for each Output in this network, i.e., 
T(y) = = {NL, NM, NS, ZN, PS, PM, PL}. The fuzzy logic 

rules for the LFLC are stated as follows: 

R t = IF x x is NL and x 2 is ZN, THEN the consequence is PM 
where i = 1, 2,..., 25. Thus, a total of 25 rule nodes are in this layer. The input to each 
rule node uf comes from one possible combination of the outputs of the Input Term 
Nodes set o} fk , where j= 1, 2,...25 and i=k= 1,2...5. Let o\ x and o| , be the possible 
inputs to the rule node Rf, however, only the smallest value of o?* becomes the input 
to the rule node Rj. Furthermore, the links (w) of this layer are unity. If there are two 
rules: 

R x =IF x x is NL and x 2 is ZN,THEN the consequence is PM (5) 
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R 2 = IF x x is NL and x 2 is ZN,THEN the consequence is PS (6) 
then the firing strengths of and R 2 are denoted as otj and (Xj , respectively. For 
example, a, is given by the following equality: 

a, = aM™(x 2 ) (7) 

where a is the fuzzy logic AND operation. Die intersection is used as the fuzzy logic 
AND operator. Dius, the AND operator is realized by the following equation: 

rmin (M «! (*i) , M q x2 {x 2 ) ) 

a, = (jcj) a M q x2 (x 2 ) =| or (8) 

M q xl ( Xl )M q 2 (x 2 ) 

where q is one of the {NL, NS, ZN, PS, PL} and i = 1, 2,...,25. Die nodes in this 
layer form a fuzzy rule base. Die connectionist inference engine is constructed by 
combining the functions of this layer and the functions of layer 4 [28], [7], [5]. Hence, 
the rule matching process is avoided [18]. Die precondition matching of fuzzy logic 
rales is accomplished by the method of linking layers. Each rale node in this layer per- 
forms the fuzzy AND operation; thus the activation function (a) for the rale node is the 
minimum of all its inputs. 

a =f (9) 

/ = mini i?, 1 (10) 

'12 p ' 

and for the network analyzed p = 2. 
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The connections of this layer are given as: a node in each of the Input Term Nodes 
set is connected to a rule node with a constraint such that no two nodes are from the 
same Input Term Nodes set are connected. 

3.5 Output Term Nodes Layer 

This a layer has two operation modes. In the first phase of the training, the nodes 
of this layer operate in an up-down transmission mode. Whereas in the second phase 
of the training, these nodes operate in a down-up transmission mode. Upon comple- 
tion of the learning process, the set M y - {M l y , M *, ..., M y } of the membership 
functions of the Output Term Nodes set are found. An Output Term Node number j 
may be excited (as a result of the learning process in which the structure of connec- 
tions between the Rule Nodes layer and the Output Term Nodes layer are established) 
by a few or none of output signals of the Rule Nodes. 

3.5.1 Down-Up Mode 

In this mode, the node performs the fuzzy logic OR operation. To illustrate this 
concept, let us assume that (as defined by the structure of connections) an Output Term 
Node is excited by the Rule Nodes number ij, i 2 ,...,/ p that is described in section 3.4. 
In our particular structure, p can be any number in the set { 1 , 2,..., 25 ) at the Rule Node 
layer; Then, the membership function associated with the Output Term Node number 
j can be defined as follows: 


M‘f J (s it ) - a it A M‘(s it ) « \ 

where i - 1, 2 


min (a, v (s it ) ) 


or 


(ID 


l V-iMy 
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where s it is the output of the rule node number i*. m y and 5 y are the mean and the 
variance of the membership y, respectively. In the case considered, we have a set of 
membership functions M y = [M], M ], ..., MJ} . Combining (8) and (11), we obtain 
the output decision: 



Mi = M l / j {s-J v M? J (sj v ... v fyjr* (s if ) (13) 

where v is the fuzzy logic OR operation which performs the UNION of a given set 
of memberships. Thus, the output decision for an Output Term Node number j can be 
written as follows: 

max (My’ j (s,-) , M‘/ j (sj M‘f j (s if ) ) 

Mi = ■ or (14) 

min(l, X fy'y J (s it ) ) 

u 1 

As a result, the nodes in this mode perform the fuzzy logic OR operation to inte- 
grate the fired rule nodes, which have the same consequence; that is, they are con- 
nected to the same Output Term Node j. This can be also expressed as 

dj = min (1, ft (15) 

where 
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f;=2Mi rl (\) (16) 

k = 1 

Each of the 25 Rule Nodes can be potentially connected to each of the seven Output 
Term Nodes as can be seen from Table 1 showing connection weights. Those weights, 
either 0 or 1, are established in the learning phase and tell us which i*’s are involved in 
equation (13) for a given node number j. 


TABLE 1. The links (weights) of the Output Term Node Layer 


Rule number 

Ml 010000000001 
M2 100001000000 
M3 000000000000 
M4 000000000000 
M5 000000000000 
M6 000010000010 
M7 0011001 1 1 100 


-25 

11110 11111111 
0000100000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
000000000000 0 



3.5.2 Up-Down Mode 

The purpose of this mode is to find the initial means and the variance of the Output 
Term Nodes. This mode is used in the first (unsupervised) learning phase. In this 
mode, the nodes in this layer, function the same as those in the second layer, except 
that only one node is used to perform a membership function for the Output Linguistic 
Node. 

3.6 Output Linguistic Nodes Layer 

This layer contains the nodes performing the UPDOWN transmission and the 
nodes performing the DOWNUP transmission. The UPDOWN transmission nodes 
are used to feed the training data into the LFLC network [18]. Their activation is 
defined as follows: 


a =f=y (17) 

The DOWNUP transmission nodes together with the links that are attached to the 
Output Linguistic Node, performs the fuzzy OR operation or in other words imple- 
ments the so called defuzzifier [18]. In this research, the fuzzy OR operation is based 
on the center of area method, which as described in [3], gives the best result. Let sj be 

the support value, i.e., a value at which the membership function, M y ( s ) , reaches the 


maximum value M y (s) j j _ j . Then, from (13) the defuzzification output is 


y = J. 


(18) 
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The following equations are used to simulate the center of area defuzzificaton method 
[32]: 


f= !>?• «? = I(m?.8?.).«f 

(19) 

/ 

(20) 


where the wfj is the link between the Output Term Node number j and the Output Lin- 
guistic Node number i, it is also equal to the product of m? • 8 ?■ = mj • 8 $ . Here m 
and 8 are the mean and the variance of the Output Term Node, respectively. In this 
particular network, we only have one Output Linguistic Node, thus j is equal to 1. 
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Chapter 4 


4.1 Learning Phase 

In LFLC there are two learning phases: the first one implements the Unsupervised 
or self-organized learning, and the second one implements the Supervised Learning 
using Backpropagation. As a result of Unsupervised Learning, the structure of the 
LFLC is established after the links (w) and the firing strengths of the fuzzy logic rule 
nodes in the network are found. Furthermore, the means (m) of each membership 
function in the network are found by using the Self-Organized learning method which 
is described by the Kohonen’s feature-map algorithm [14], [13]. 

4.2 Unsupervised Learning 

The Kohonen feature-map algorithm is a self organizing method that gathers the 
input data into a cluster. In his algorithm, a set of weights is initially generated for 
each node in his network. When data is present in the network, the distances from the 
data to all nodes are calculated. Let w i=i,j=i denote the weights of the shortest dis- 
tance from the data i to the node j, then the node number j is selected to be the output 
node of the network. Furthermore, the weights (w ,_j J=J ) and its neighbors are 
updated according to the following equation: 

Wij (r + 1) = Wij (t) + r\ ( t ) (Xi (0 - wu ( t ) ) (21) 

where j belongs to the nearest neighbor, and 0<i <N- 1 . The term rj (t) is a learn- 
ing rate (Ocq (t) < 1 ) which decreases in time. 
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After the mean of the Fuzzy Logic membership function is found, the correspond- 
ing variance (5) can be found as follows: 

( 22 ) 

where x is the overlap parameter. This parameter is used to control the level of over- 
lapping between memberships in the same cluster. The range of these parameters 
depends on the range of the input data. However, in our network, this parameter is 
chosen to be 1.5. 

The m c[osest is the nearest neighbor of the current mean value. Once the centers 
and the variances are found, the input signal and the teaching signal reach the output 
points at the output term nodes and the input term nodes. Next, the output of the input 
term node at layer two are transmitted to layer three through the initial architecture of 
the layer three’s links. Based on the firing strengths of the rule nodes (output of the 
rule nodes) and the output of the output term nodes, the correct consequences of the 
links of each fuzzy logic rule node are found. Initially, the links (w) are fully con- 
nected. However, the competitive learning algorithm is used to update the links 
(weights) for each training data set This algorithm is described by the following 
equation: 

Wij (t) = of (- + of) (23) 

where of, of , are the outputs of the output term node and the input term node, respec- 
tively. w i} denotes the weight of the link between the /-th rale node and the y'-th output 
term node. A dot over w 0 denotes the next value of the w {j . 


g _ \ m i ~ m cloteti 
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After competitive learning, the weights of the links at layer four represent the fir- 
ing strength of the rule node which is transmitting the consequence of the rule node to 
the term node of an output linguistic node. Furthermore, the links are chosen such that 
at most one link is selected and the others are eliminated. As a result, only one term 
node in the output term node set becomes one of the consequences of the fuzzy logic 
rule. The supervised learning takes place after the fuzzy logic rule is established in the 
unsupervised learning phase. 

4.3 Supervised Learning 

In the supervised learning phase, the backpropagation method is used to fine tune 
the parameters of the LFLC which were described above. The objective of the second 
phase of learning is to minimize the error function: 

E = 3 (y(r) -?(*))* (24) 

The error signals for layer five, layer four, and layer three, are given as: 


a 5 (r) = y(t) -y(t) 


erf (t) 


m t°i ( £5 ;“.) - (EmfiiUi) 5, 
W ‘ (L5 u) 2 


(25) 

(26) 


cf = Xo? (27) 

k 

where y(t),y(t) are the desired output and the current output of the LFLC, respec- 
tively. 

Backpropagation has a forward phase and a backward phase. The forward phase 
of the backpropagation requires that the data be presented to the network at the first 
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layer for training. Next, the outputs of the input term nodes are calculated, and then, 
transmitted to one of the output term nodes, via the firing strengths of the rule nodes. 
Each fuzzy logic rule was structured to be excited by the smallest output value of the 
initially defined input term nodes. Finally, the output of the LFLC is calculated using 
equation (18). After the output is found for a pair of the training data, the error signal 
is calculated using equation (24) and is then propagated to all of the previous layers of 
the LFLC. Concurrently the error signal is being used to update the set of means (m) 
and the set of variances (8 ) of each layer in the LFLC.The means and the variances at 
layer five are updated (fine tuned) by equations (27) and (28), respectively. 


m? (t + 1) = mf (f) + T)-o 5 - 

mu (25 ,«,) - (Sm.-S.u,) u ; 


8f(r+ 1) = Sf(t) + n*c 5 


(Z5 iUi ) : 


(28) 

(29) 


where the learning rate (Tj ) is a function decreasing in value as the time progresses. 
Here the symbol u is the input value to the current node i and 5, m is the variance, and 
the mean of the current membership i, respectively. Layer four does not contain any 
parameters. Thus, only the error signal (c 4 ) of this layer needs to be calculated, and 
it is then propagated to the layer three. The equation of the error signal (cr 4 ) is 
obtained from equation (25). As shows in Fig. 1., layer three also does not contain any 
parameters that have to be updated, hence only the error signal (a 3 ) is calculated and 
propagated to layer two. In the layer two, the mean (m) and the variance (5) are 
updated according to the following: 
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The general learning rule in Neural Network is 

Awoc(-M) (30) 

which can be written as 


W ( I+ i) = 

where rj is the learning rate and E is the error defined by eq. (24) and 


(31) 


BE BE d (net -input) _ BE df 

- B (net -input) Bw ~ Bf Bw 

where 


(32) 


net-input = f(uf, u%, ...,u*;w},w£, ...,w*) . 

here, the superscript k indicates the layer number and subscript p indicates the index of 
the input nodes. 


BE 

Bw 


BE Ba Bf 
Ba ’ Bf Bw 


(33) 


BE 

To show the learning rule, we will derive the computation of for this layer. 


Using (32), and (4), the adaptive rule for m,* is derived as follows: 


BE _ BE Ba; Bf_ _ BE , 2 (u,-m u ) 
B(m u ) ~ Ba, Bf Bm v Ba , o?- 


(34) 
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where i is the number of the node in the Input Term Node layer, and j is the number of 
the node in the Rule Node layer. 

We have 




3E 3 (net- input) * 

k 3 ( net - input ) * 3a,- 


(35) 


for this particular layer; 


3 E 


3 E 


o* 


(36) 


3 (net - input) * djf 
where o% is the error signal for rule node at layer three, and k= 1, 2,...,25. Also from 
(9), and (10) 

d ( net l,i/ «|= min( inputs of rule node k), otherwise 0. 


Hence, 



(37) 


where the summation is performed over the rule nodes that a,- feeds into. Here a,- is the 
output of the Input Term Nodes number i and 


f c l , ifai is minimum in kth rule node ’s input 

(38) 

0 , otherwise 

Thus, the adaptive rule of m,y for this layer is 
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m y (t+ 1 ) = m,j (t ) ~ ^ 

In the same manner, using (32), (4) and (37), (35),the adaptive rule of 5 0 is 




■ J 




i 
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Chapter 5 


5.1 Symbolic Representation and Performance Evaluation of the LFLC 

In the LFLC, the intelligent control decision of the controller is determined by 
both the links (or weights) of the Learning Fuzzy Logic Controller and the variables of 
the Fuzzy Logic membership functions. The symbolic representation of the network 
is described in two levels: the Low Level, and the High Level. 

5.2 Low Level Symbolic Representation of the LFLC 

The low level symbolic representation of the network consists of the variables of 
the fuzzy membership function and the links in each layer of the network. We can 
interpret the symbolic representation of this level as the IF part of the LFLC. For 
example, the mean is the center of a set of data that belongs to a cluster. The cluster 
method that is used in this network is Self-Organization and the variances of the fuzzy 
logic membership function serve as the radius of the cluster. Data that is out of the 
range of the radius will have no effect on the node in either the Input Term Nodes layer 
or the Output Term Node layer. Thus, the combination of the mean and the variance of 
the Fuzzy Logic membership function as an activation function of a node in the neural 
network enables the network to perform decision making. In an extreme case, each 
node in the network represents Off or On in linguistic terms for the LFLC. The Fuzzy 
Logic membership functions following an unsupervised learning process are shown in 
Fig. 3 and Fig. 4. Each Gaussian membership function has a mean (center) and a vari- 
ance (width). If data are presented to these Fuzzy Logic membership functions, only 
one membership function will have the highest activation value. 
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Please note that the membership function 1-2 is not shown in Fig. 3 and Fig. 5. 
This is due to the fact that the mean and variance values of this membership function 
did not contribute to the control decision making (the mean value and variance did not 
change their values) for this particular node. Hence, it is not shown in the Fig. 3 and 
Fig. 5. 



Input Data (xl) 


■ Membership 1-1 + Membership 1-3 * Membership 1 -4 □ Membership 1 -5 


Fig.3. Initial Input Term Node membership for xj 
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Input Data (x2) 

-m- Membership 2-1 — Membership 2-2 -**- M ember shp 2-3 
hb- Membership 2-4 -M- Membership 2-5 

Fig.4. Initial Input Term Node membership for X 2 

After supervised learning, most of the mean and variance values of the Fuzzy 
Logic membership functions are changed from their previous values (before the super- 
vised learning). This is shown in Fig. 5 and Fig. 6. In Fig. 5, the mean and variance 
values of the Fuzzy Logic membership functions are different than those of the mem- 
bership functions shown in Fig. 3. In the same manner, the mean and variance of 
membership functions 2-1, 2-2, 2-3, 2-4 and 2-5 in Fig. 6 are different than those of 
the membership functions shown in Fig. 4. This is due to the fact that in the super- 
vised learning phase, the error between the calculated value and the desired value is 
calculated and propagated to all layers in the network. As a consequence, the mean 
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and the variance are updated according to the learning rule (30). After supervised 
learning, the mean and the variance represent the best control decision for a specific 
node. 



■ Membership 1-1 ♦ Membership 1-3 w Membership 1-4 a Membership 1 -6 


Fig.5. Final Input Term Node membership for Xj 


; 


Also, the variances of the Fuzzy Logic membership functions 1-3 and 1-4 in Fig. 3 
changed their values from large value to small that is shown in Fig. 5. This can be 
seen better when the figure is enlarged. This indicate that only limited range of data 
will be in the neighbor hood of the means of these two membership functions. Other- 
wise, the variances of these two membership functions would increase in value, like 
those of the membership functions 1-1, and 1-5, to allow more data excites the mem- 
bership functions that gives higher activation value. Also, the means of these two 
Fuzzy Logic membership functions move closer to each other. The changing of direc- 
tion indicate the input data that belongs to these two clusters is somewhat separate 
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from other input data. Furthermore, the variances of the membership functions 1-1, 
and 1-5 are increased in value. This increase in value indicate that majority of the data 
belongs to these two clusters, and the center of these two clusters are the mean values 
of those two Fuzzy Logic membership functions 1-1, 1-5. 



InpLt Data {x2) 


~m- Membership 2-1 Merriwshp 2-2 -*f- Membership 2-3 
-e- Membership 2-4 -w- Menfcershp 2-5 

Fig.6. Final Input Term Node membership for x 2 
In general, a change in the mean indicates that most of the data is close to the new 
value. A large value of variance indicates that a large quantity of data is located at that 
mean, and conversely, smaller values of the variance indicates that smaller quantities 
of the data are at that mean value. The convergence of the mean and the variance are 
shown in Fig. 7, 8, 9, and 10. 
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Standard Deviation 
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• Mean 2-1 Mean 2-2 ** Mean 2-3 
■ Mean 2-4 — — Mean 2-5 


Fig.9. Convergence of the means of the Input Term Node memberships for input x 2 



Nurrtxsr of Iterations 


— Deviation 2-1 Deviation 2-2 m Deviation 2-3 
Deviation 2-4 Deviation 2-6 

Fig. 10. Convergence of the standard deviations of the Input Teim Node memberships for input x 2 
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In a Neural Network, the links represent the firing strengths which are used to con- 
nect the nodes of one layer to the other layers. However, three types of links are used 
in LFLCS. The first type of link is used to transmit information from one node to 
another node. For example, the links of the Input Linguistic Node layer are used as a 
connection to transmit information from this layer to a node at the Input Term Node 
layer. The connection of the network is shown in Fig. 2. In the same manner, links at 
the Output Term Nodes are used to transmit signals from the nodes in this layer to the 
nodes at of Output Linguistic layer, and the links of these layers are fully connected 
with weights equal to 1. 

In this particular network, the second type of link is constructed in such a way that 
for each node in the Rule Nodes layer there are only two connections from the Input 
Term Nodes layer. These links can then be interpreted as the IF part of the Fuzzy 
Logic rule. The predefined structure of these links is explained in detail in section 3.4. 

The third type of link used is at the Output Term Node layer. These links are ini- 
tially fully connected, but the weights are modified by the learning rule (eq. 23) in the 
first (unsupervised) phase of the learning process. These links can be interpreted as 
the THEN (consequence) part of the Fuzzy Logic rule. As can be seen in Table 1 , a 
rule node can be connected to only one node at the Output Term Nodes layer. These 
links represent the firing strength for a node in that network. Furthermore, Table 1 
shows that nodes number 3, 4 and 5 of the Output Term Node layer are eliminated. 
This is because none of the rule nodes connects to these Output Term nodes. Thus, 
any rule node that contains one of the three nodes from the Output Term Node, men- 
tioned above as its consequence node, is automatically eliminated. 
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5.3 High Level Representation of the LFLC 

In Fuzzy Logic symbolic representation of control strategy, the use of a linguistic 
term offers advantages over the conventional approach to specifying the control algo- 
rithm as an equation, especially, in ill-structured situations. The concept involves 
using a linguistic rule to describe the operation of the process from a human point of 
view and to capture the essential knowledge of the operation of that process, which the 
operator has presumably acquired through direct experience with and actual operating 
process. It follows that this knowledge can be used as the best rule set which the oper- 
ator can obtain for the control action in linguistic terms. 

In this thesis, the linguistic term of the control action is learned through the train- 
ing of the network, and the teaching pattern is used for operator’s knowledge. This 
knowledge is stored in the network by first using the Self-Organizing learning method 
and then using the Backpropagation method of learning in order to determine the 
Fuzzy Logic rule nodes. 

For example, the linguistic term for the first node in the Input Term Nodes layer is 
interpreted as Negative Large (NL). The second node is interpreted as Negative Small 
(NS), the third node is interpreted as Zero ( ZN ), the fourth node is interpreted as Posi- 
tive Small (PS), and the fifth node as Positive Large (PL). 

The linguistic term for the first node of the Output Term Nodes is interpreted as the 
Negative Large (NL) node. The second node is interpreted as Negative Medium (NM), 
the third nodes can be interpreted as Negative Small (NS), the fourth node is inter- 
preted as Zero (ZN), the fifth node is interpreted as Positive Small (PS), the sixth node 
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Fig. 1 1 . Initial Output Term Node membership 
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5.4 Performance evaluation 

In this section, performance of the LFLC is evaluated by an example using the 
LFLC in a closed loop control system. First, a conventional PI controller is designed 
to control the shaft position of a D.C. motor. The transfer function of the D.C. plant 
G p (z) is 


1 79,45x1 0 -6 (z + 1) 

G ' u,= (I^Tp 


(40) 


The transfer function of a position sensor H(z) which is used in the feedback path of 
the closed loop system is 


H(z) = 


99. 3z - 95.3 
z + 1 


A unit step signal is used as the input signal to the system. 


(41) 



— — Control Signal Error Signal Feedback Signal 


Fig.13. Training Data 
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The control system is implemented by simulation to obtain training data set as 
shown in Fig. 13. This data set is then used in the off-line training of the LFLC to 
obtain the fuzzy logic controller, which will give approximately the same performance 
as a PI controller. 

A closed loop control system is also simulated to illustrate the capability of the 
LFLC. This closed loop control system is shown in Fig. 14. In the diagram, the error 
signal and the feedback signal are normalized before being fed into the LFLC. This 
normalization process is a necessary step before any data is presented to the LFLC. 
This follows since during the learning process the training data is also normalized by 
equation (1). However, it should be observed that the control signal generated from 
the LFLC in Fig. 14 is going through an inverted normalization processes in order to 
retrieve the actual control signal. This is required since the teaching pattern is normal- 
ized by equation (2) during the training process. Equation (43) is used in the inverse 
normalization process. 



Fig.14. LFLCS with PI Fuzzy Logic Controller 


( 42 ) 
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y actual - (yiFic-0-5) • (2xmax(\maxy\,\miny,\)) (43) 

where y^c is the output signal from the LFLC, and y actual is the actual control sig- 
nal. 

The output of the computer simulation of the LFLC with a PI Fuzzy Logic control- 
ler is shown in Fig. 15. The result of the simulation shows that the LFLC is able to 
control the D. C. motor at least as good as a conventional PI controller or even better. 
The LFLC shows that it has better rise time performance, better settling time and a 
smaller steady state error. 



— LFLNC MaUab Simulation 


Fig. 15. Step response of the system with classical and fuzzy logic PI controller 
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Chapter 6 


6.1 Summary 

In this thesis, the Learning Fuzzy Logic Controller is being developed and replaces 
the conventional PI controller that is used in a closed loop position control system. 
The result of two simulated control systems is given in Fig. 15. We find that the rise 
time, settling time and steady state error of the system with the LFLC are superior over 
those of a system with a conventional PI controller. 

The Learning Fuzzy Logic Controller developed in this thesis shows that by com- 
bining the Neural Network learning concept and the Fuzzy Logic rule base, it is possi- 
ble to eliminate the need for an accurate model of the system. 

However, during the research of the LFLCS we discovered that there are some 
limitations of this approach. This is due to the need of the teaching pattern necessary 
to train the network. In most cases, this teaching pattern is the control signal to the 
plant This signal can be either generated by computer simulation as we did in this 
thesis or using an actual control signal of a real system. However, in control applica- 
tion, this signal is not always obtainable. In order to overcome this problem in control 
design, Reinforcement learning approach for control design should be considered. 

The reinforcement learning utilizes both the knowledge of the system at hand and 
error prediction scheme to modify the parameters of the controller. 
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APPENDIX B 
Simulation Code 


#include <stdio.h> 

#include <stdlib.h> 

#include “c:\borlandc\file\Xhesis.h” 
#include“c:Nborlandc\file\globalex.h” 


ji'l 


int 

mainO 

{ 


unsigned int input_ling_no, 
output_ling_no, 
i, 
j; 


if (!open_files()) 

{ 

printf(“CAN NOT OPEN FILENn”); 

> 


if (!read_input()) 

{ 

printf(“CAN NOT READ INPUT FILEmT); 

} 


if (!unsupervised_leaming()) 

{ 

printf(‘TNITIAL PHASE1 FALSEXn”); 
return (FALSE); 

> 


if 

K-V.V 


if (!supervise_leaming()) 

{ 

printf(“BACKPROP PHASE1 FALSExn”); 
return (FALSE); 

} 

printf(“end of supervise leaming\n”); 
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getchar(); 

* A 

if (!close_loop_control()) 

{ 

printf(“CLOSE CONTROL LOOP FALSER”); 
return (FALSE); 

> 

return (TRUE); 
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#include <stdio.h> 

#include “c;Mx>rlandc\fileVhesis.h” 

#include “c:'iiorlandc\file\globalex.h” 

in t 

open_files(void) 

{ 

if ((Error_fptr - fopen(“error.txt’\ “r”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN ERROR FELENn”); 
return (FALSE); 

> 

if ((Fback_fptr •* fopen(“fback.txt”, “r”)) — NULL) 

{ 

fprintf(stdeir, “CAN NOT OPEN FBACK FILE'n”); 
return (FALSE); 

> 

if ((Consigna_fptr - fopen(“consigna.txt”, “r”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN CONSIGNA FTLEmT); 
return (FALSE); 

> 


if ((Output_fptr - fopen(“outputtxt”, “w”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN OUTPUT FILE'n”); 
return (FALSE); 

> 

if ((W_I4_fptr - fopen(“weigth4.txt”, “w”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN WEIGTH4 FELEvn”); 
return (FALSE); 

> 

if ((Input_term_Deviation_fptr - fopen(“input_D.txt”, “w”)) — NULL) 

( 

fprintf(stderr, “CAN NOT OPEN INPUT DEVIATION FILEn”); 
return (FALSE); 

> 

if ((Fina]_input_term_deviation_fptr - fopen(“fi JnD.txt”, “w”)) — NULL) 
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{ 

fprintf(sjtderr, “CAN NOT OPEN FINAL INPUT DEVIATION FELEmT); 
return (FALSE); 

> 

if ((Output_term_Deviation_fptr - fopen(“output_D.txt”, “w”)) — NULL) 

{ 

fprintf(stdeiT, “CAN NOT OPEN OUTPUT DEVIATION FILEmi”); 
return (FALSE); 

> 

if ((Final_output_term_deviation_fptr - fopen(“fi_outD.txt”, “w”)) NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN FINAL OUTPUT DEVIATION FILE\n”); 
return (FALSE); 

> 


if ((Mean_output_fptr - fopen(“mean_out.txt”, “w”)) — NULL) 

{ 

fprintf(stderr, "CAN NOT OPEN MEAN OUTPUT FILEXn”); 
return (FALSE); 

> 

if ( (Final_mean_output_fptr - fopen(“fin_mout.txt”, “w”)) -- NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN FINAL MEAN OUTPUT FELEvn”); 
return (FALSE); 

> 


if ((Mean_input_fptr - fopen(“mean_in.txt”, “w”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN MEAN.IN FILE\n”); 
return (FALSE); 

> 

if ((Final_mean_input_fptr - fopen(“fin_min.txt”, “w”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN FINAL MEAN INPUT FELESn”); 
return (FALSE); 

> 

printf(“OPEN SUCCESSFULV’); 
return (TRUE); 

> 
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#include <stdio.h> 

#include <stdlib.h> 

#include “c:\borlandc\file\thesis.h” 
#include‘‘c:\borlandc\file\globalex.h” 


* function init_phasel() 

* 

* This function excute all initial procedure for getting the center of the 

* membership and the deviation of input and output term nodes.So after this 

* function the input signal can be pass through the membership on layer3 and 

* layer4. 

* 

♦Input parameter, inputxl, 

*inputx2. 

* 

* Return value: TRUE if each of the center of membership and the deviation 
♦are calculated. 

♦FALSE otherwise. 

* 

* Programmer 

* Leung Kam Lung 9/30/93MS. Thesis. 

* 

*/ 

int 

unsupervisedlearning(void) 

{ 

inti; 

float time, 
time2; 

unsigned int input_lingui_no, 
output_lingui_no; 

time -0.001; 
time2- 0.001; 
initials_all(); 
input_lingui_no - 0; 
output_lingui_no - 0; 

while ((time > 0.000001) II (time2 > 0.000001)) 

{ 

for (i - 0; i < MAX_INPUT_INDEX; i++) 

{ 

update_input_12_w(Inputx[input_lingui_no][i], time, input_lingui_.no); 
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update_input_12_w(Inputx[input_lingui_no + l][i], time, input_lingui_no + 1); 
upda(e_output_w(Output[outputJingui_no][i], time2, output_lingui_no); 

> 

time - time * 0.995; 
time2 *» time2 * 0.995; 


> 


if (!save_mean_input()) 

{ 

printf(“Mean_input FALSER”); 

} 

if (!save_mean_output()) 

< 

printf(“Mean_output FALSEmi”); 

> 


/* 

* up to now only the center(mean) of the input and output membership are 

* calculated 
*/ 

for (input_lingui_no - 0; input_lingui_no < INPUT_LINGUI_NO; input_lingui_no++) 

< 


if (find_input_deviation(input_lingui_no)) 

{ 

printf(“ Deviation of each Membership is find for INPUT_TERM_NODE %<J\n”, 
input_lingui_no); 

> 

else 

{ 

printf(“ NOT Deviation of each Membership is find for INPUT_TERM_NODE %d\n”, 
input_lingui_no); 

return (FALSE); 

> 


} /* End for input_lingio_no */ 

if (!save_input_deviation()) 

{ 

printf(“CAN NOT SAVE INPUT DEVIATION^”); 

> 
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for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


if (find_output_deviation(output_lingui_no)) 

{ 

printf(“ Deviation of each Membership is find for OUTPUT_TERM_NODE %d\n”, 
output_lingui_no); 

> 

else 

{ 

printf(“ NOT Deviation of each Membership is find for OUTPUT_TERM_NODE %d\n”, 
output_lingui_no); 

return (FALSE); 

} 


> /* End for output_lingio_no */ 

if (!save_output_deviation()) 

{ 

printf(“CAN NOT SAVE OUTPUT DEVIATIONS”); 

> 

/* 

* now all deviation of each membership for input term nodes and output 

* term nodes are calculated 
*/ 


if (!update_w_14()) 

{ 

printf(“udate_w_14 FALSENn”); 
return (FALSE); 

> 

/* printf(“ after updata w 4\n”); */ 

if (!find_max w4()) 

{ 

printf(“can not find max weight 4 FALSEs”); 
return (FALSE); 

> 

/* printf(“ after find_max_w4 \n”); */ 
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if (!save_w4{)) 

< 

printf(“can not save weight 4\n”); 
return (FALSE); 

> 

/* printf(“ after save w 4\n”); */ 


return (TRUE); 



#include <stdio.h> 

#include <st$lib.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

int 

supervise_learning(void) 

{ 

float Irate; 
int index, 
retv5, 
retv4; 

Irate - 0.045; 

/* 

* forward phase of the learning for all inputs. At this point all 

* Errors will be calculated and then updata procedures will then 

* execute. Error will be sumed for smoth result 
*/ 

Cal_error « 1.0; 

while ((Irate > MINJLRATE) && (Cal.error > 0.025)) 

{ 


Cal_error - MIN_ERROR; 

Max_indexs - 0; 

for (index - 0; index < MAX_INPUT_INDEX; index++) 

{ 


/* printf(“NOW IN THE BACKPRO FORWARD \n”); */ 
if (!backpro_forward(index)) 

{ 

printf(“ Can not find the output from backpro forwardsn”); 
return (FALSE); 

} 

/* printf(“NOW IN THE BACKPROV’); */ 

/* getcharO; */ 

if (!error_backpro_layer5(lrate, index)) 

{ 

printf(“ Can not find the error_backpro_layer5_ptiNn”); 
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return (FALSE); 

> 


/* printf(“\n after LAYER 5\nW’); */ 

if (!error_backpro_layer4(Sum_width_act, Sum_mean_width_act)) 

{ 

printf(“ Can not find the error_backpro_layer4jpti\n”); 
return (FALSE); 

> 

/* printf(‘V after LAYER 4\nV’); */ 

if (!error_backpro_layer3()) 

{ 

printf(“ Can not find the error_backpro_layer3_ptf\n”) ; 
return (FALSE); 

> 

/* printf(“\n after LAYER 3\n\n”); */ 

if (!error_backpro_layer2(lrate, index)) 

{ 

printf(“ Can not find the error_backpro_layer2_ptf\n”); 
return (FALSE); 

> 

/* printfCV after LAYER 2\n\n”); */ 

/* printf(“\nNOW UPDATE WEIGHTW’); */ 

if (!update_w(index)) 

{ 

printf(“Error update the weighin’’); 
return (FALSE); 

> 

/* getchar(); */ 

) /* End of for input_index */ 

Irate *- 0.998; 

printf(“ actual output %f Output[0] [Max_indexs]); 
printf(“ calculated output %fV’, 

Output_down_up[0] [Max_indexs]); 

printf(“leam rate % 1 . 1 Of Irate); 
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> 


printf(“Cal_error %1.10f\n”, Cal_error); 

' ,<i 

f* End of the while Irate loop */ 

save_final_mean_input(); 

save_final_mean_output(); 

save_final_input_deviation(); 

save_final_output_deviation(); 


return (TRUE); 

> 
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#include <stdio.h> 

#include <sl<Jlib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


int 

get_downup_act4(unsigned int input_index) 

{ 

unsigned int output_lingui_no; 
unsigned int output_term_membership; 
unsigned int rule_index; 
float act4; 

for (output_lingui_no « 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 
output_terxn_membership++) 

{ 


act4 = 0.0; 


for (rule_index «* 0; rule_index < RULE_NODE; rule_index++) 

{ 

act4 +- Act_13[rule_index] * 

W_14[output_lingui_no] [output_term_membership] [ruie_index] ; 

> /* End rule_index */ 


if (act4 >= 1) 

{ 

Act_downup_14[output_iingui no][output_term_membership] - 1.0; 

> 

else 

{ 


Act_downup_14[output_lingui_no][output_term_membership] - act4; 
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> . 

* it 

if (act4 > MAXFLOAT) 

{ 

printf(“activatio value of level 4 is bigger that max floatSn”) 
return (FALSE); 

> 

} /* End for output_term_membership *1 


> /* End for output_lingui_no */ 

return (TRUE); 
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#include <stdio.h> 

#include <st$ilib.h> 

#include <math.h> 

#include “c:\borlandc\file\lhesis.h” 
#include “c:\borlandc\file\globalex.h” 


int 

get_downup_act4(void) 

{ 

unsigned int output_lingui_no; 
unsigned int output_term_membership; 
unsigned int rule _index; 
float act4; 

for (output_lingui_no ** 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 

output Jerm_membership++) * 

{ 


act4 - 0.0; 

for (rule_index - 0; rule_index < RULE_NODE; rule_index++) 

{ 


act4 +- Act_13trule_index] * 

W_14[output_lingui_no] [output_term_membership] [rulejndex] ; 
} /* End rule_index */ 


if (act4>“ 1) 

{ 

Act_downup_14[output_lingui_no][output_term_membership] - 1.0; 

} 

else 

{ 


Act_downup_14[output_lingui_no][output_term_merabership] - act4; 
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} 


if (act4 > MAXFLOAT) 

{ 

printf(“activatio value of level 4 is bigger that max fioatNn”); 
return (FALSE); 

} 

} /* End for output_term_membership */ 


> /* End for output_lingui_no */ 

return (TRUE); 

> 
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#include <stdio.h> 

#include <stfllib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

/************************** ************** ** *********** *** ****** * *** ************* 

* function backpro() 

* 

* This function calculates the network output in downup pass of the network. 

* It call the following functions. 

* 

* get_downup_act4; 

*get_act3; 

*get_act2; 

* 

* 

* Input parameter, index 

* 

* Return value: int type. 

★ 

*programmer:Leung Kam LungMS. Thesis9/28/93 
* 

*/ 


int 

backpro_forward(unsigned int index) 

{ 

int retval2, 
retval3, 
retval4; 

unsigned int input_lingui_no, 
input_lingui_nol, 
input_lingui_no2, 
output_lingui_no; 

float outl, 
f, 

a, 

temp_error, 

abs_cal_error, 
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unsigned int output_term_membership, 
membership; 

input_lingui_nol - 0; 
input_lingui_no2 « 1; 
output_lingui_.no - 0; 


/* printf(“ index at backpro %d\n”, index); */ 

for (input_lingui_no - 0; input_lingui_no < INPUT_LINGUI_NO; input_lingui_no++) 

{ 

retval2 = get_act2(input_lingui_no, index); 

if (retval2 — FALSE) 

{ 

printf(“reture value of act_function_ptr2 is > or < max, min doubleW’); 
return (FALSE); 

> /* End of if retval2 */ 


} /* End of for input_lingui_no */ 

retvaB - get_act3(input_lingui_nol, input_lingui_no2); 

if (retval3 — FALSE) 

{ 

printf(“reture value of act_function_ptr3 is > or < max,min doubleNn”); 
return (FALSE); 

> /* End of if retval3 */ 


retval4 - get_downup_act4(); 

if (retval4 — FALSE) 

{ 

printf(“reture value of act_function_downup_ptr4 is > or < max, min doubleW’); 
return (FALSE); 

> /* End of if retval *! 
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/* 

* The next for loop calculates the network’s output of the supervise 

* learning 
*/ 

for (output_lingui_no = 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 

Output_down_up[output_lingui_no] [index] - 0.0; 
f-f2-0.0; 


for (output_term_membership = 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 
output_term_membership++) 

{ 


f +« (Mean_output[output_lingui_no][output_term_membership] * 
Output_term_Deviation[output_lingui_no][output_term_membership] * 
Act_downup_14[output_lingui_no][output_tenn_membership]); 

f2 +- (Output_term_Deviation[output_lingui_no][output_term_membership] * 
Act_downup_14[output_lingui_no][output_term_membership]); 

> /* End of output_term_membership */ 

outl - Output_down_up[output_lingui_no] [index] - f / f2; 

Error_array[output_lingui_no] - temp_error - (float) (Output[output_lingui_no] [index] - 
Output_down_up[output_lingui_no][index]); 

/* printf(“ actual output %f“, Output[0][index]); *1 

/* 

* printf(“ calculated output %f, Error %f\n”, outl, 

* EiTor_array[output_lingui_no]); 

*1 

/* getchar(); */ 

/* abs_cal_error - (float) pow(fabs(temp_err«r) t 2); */ 

abs_cal_error - (float) fabs(temp_eiror / Output[output_lingui_no] [index]); 

if (abs_cal_error > Cal_error) 

{ 

Cal_error - abs_cal_error, 
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Max_indexs = index; 

* ft 

> 


} /* End of output_lingui_no */ 


return (TRUE); 

} 


63 



#include <stdio.h> 
#include <stfllib.h> 
#include <math.h> 


#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


int 

change_w4(void) 

{ 

float weight4; 

unsigned int output_lingui_.no; 
unsigned int output_term_membership; 
unsigned int rulejndex; 

for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_li: 

{ 

for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 
output_term_membership++) 

{ 

for (rule_index - 0; rulejndex < RULE_NODE; rule_index++) 

{ 


/* 

* printf(“output_term_membership %d, rule 

* %d’\output_term_;membership, rule_index); 

*/ 

/* 

* printf(“ weight4 %f 

* “,W_14[output_lingui_no][output_term_membership][rule_index 
*]); 

*/ 

weight4 - WJ4[output_lingui_no][output_term_membership][rule_index] « 
Act_14[output_lingui_no][output_term_membership] * 

(Act_13 [rulejndex] - 

WJ4[output_lingui_no][output_terai_membership][rule_index]); 

/* printf(“weight4 %f\n”,weight4); */ 


i_no++) 
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if (weigh t4 > MAXFLOAT) 

{ 


printf(“weight4 value of output_lingui_no %d output term membership %d rule_index &d 
is > or < expected value\n”, 

output _lingui_no, output_term_membership, rule Jndex); 
return (FALSE); 

> 


} /* End rule_index */ 

} /* End for output_term_membership */ 

> /* End for output_lingui_no */ 

return (TRUE); 

} /* End of function change_w4 */ 
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#include <stdio.h> 

#include <stdlib.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globaIex.h” 


int 

dose_loop_control(void) 

{ 


int index, 
time; 

float inputxl, 
inputx2, 

Fback; 
float y, 

yi. 

y2; 

float ul, 
u2; 

float fb, 
fbl; 

float close_inputl, 
close_input2; 
float inv_control_signaI; 
float plant_output, 
reference; 


yl - 0.0; 
yl ■» 0.0; 
ul - 0.0; 
u2 - 0.0; 
fb “ 0.0; 
fbl - 0.0; 


reference - 1.0; 

plant_output - 0.0; 

close_inputl - reference - fb; 
close_input2 - fb; 

Close_loop_input[0] - (close_inputl * Input_slopel); 
Close_loop_input[ 1 ] - (close_input2 * Input_slope2); 
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printf(“IN CONTROL LOOP Nn”); 

* .1 

for (time - 2; time < MAX_TTME; time++) 

{ 


printf(“ Error %8f\n”, close_inputl); 
printf(“ Fback %8f \n”, close_input2); 
printf(“ Input_xO %8f \n”, Close_loop_input[0]); 
printf(“ Input_xl %8f\n”, Close_loop_input[l]); 


if (!close_loop„backpro_forward()) 

{ 

printf(“ Can not find the output from close loop backpro forwardSn”); 
return (FALSE); 

} 


inv_control_signal - ((Control_signal - 0.5) / Output_slopel); 

printf(“Teach %1.6f\n”, inv_control_signal); 

plant_output - 2 * yl - y2 + 179.45e-6 * (ul + u2); 

printf(“output of the plant %f\n”, plant_output); 

fb = 99.3 * plant_output - 95.3 * yl - fbl; 
printf(“ fback %f \n”, fb); 

close_inputl - reference - fb; 
close_inpUt2 - fb; 

/* 

* printf(“Close inputl bigger %f\n”,close_inputl*Input_slopel); 

* printf(“Close input2 bigger %f\n”,close_input2*Input_slope2); 
*/ 

if ((close_inputl * Input_slopel) > Max_inputl) 

{ 

Close_loop_input[0] - Maxjnputl; 

> 

else 

{ 
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if ((close_inputl * Input_slopel) < Minjnputl) 

{ 

Close_loop_input[0] -Min_inputl; 

} 

else 

{ 

Close_loop_input[0] = (closejnputl * Input_slopel); 

> 

} 

if ((close_input2 * Input_slope2) > Max_input2) 

{ 


printf(“Max_mput2 %f \n”, Max_input2); 
Close_loop_input[l] - Max_input2; 


> 

else 

{ 

if ((close_input2 * Input_slope2) < Min_input2) 

{ 

printf(“Min_input2 %f\n”, Min_input2); 

Close_loop_input[l] - Min_input2; 

} 

else 

{ 

CloseJoop_input[l] - (close_input2 * Input_slope2); 

> 

> 

u2-ul; 

ul = inv_control_signal; 
y2-yl; 

yl - plant_output; 
fbl -fb; 

/* 

* printf(“ u2 %8f V,u2); printf(“ ul %8f\n”,ul); printf(“ y2 %8f 

* Nn”,y2); printf(“ yl %8fV’,yl); printf(“ fbl %8f\n”,fbl); 

*/ 

getchar(); 
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fclose(Output_fptr); 

fclose(Max_data_fptr); 

return (TRUE); 

> 



#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

/***************************************************** ************* ************* 

* function close_loop_backpro_forward () 

* 

* This function.calculates the network output in downup pass of the network. 

* It call the following functions. 

* 

* close_get_downup_act4; 

*close_get_act3; 

*close_get_act2; 

* 

* 

* Input parameter index 

* 

* Return value: int type. 

* 

*programmer:Leung Kam LungMS. Thesis9/28/93 
* 

*/ 


int 

close_loop_backpro_forward(void) 

{ 

int retval2, 
retval3, 
retval4; 

unsigned int input_lingui_no, 
input_lingui_nol, 
input_lingui_no2, 
output_lingui_no; 

float outl, 
f, 

12 , 

temp_error, 

abs_cal_error. 
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unsigned int output_term_membership, 
membership; 



PI 

iiSil 


input_lingui_nol - 0 
input_lingui_no2 = 1 
output_lingui_no - 0 


for (input_lingui_.no - 0; input_lingui_no < INPUT_LINGUI_NO; input_lingui_no++) 

{ 


retval2 - close_loop_get_act2(input_lingui_no, Close_loop_input[input_lingui_no]); 

if (retval2 = FALSE) 

{ 

printf(“reture value of close_loop_get_act2is > or < max,min doubleW’); 
return (FALSE); 

> /* End of if retval2 */ 

} /* End of for input_lingui_no */ 

retval3 - close_loop_get_act3(input_lingui_nol, input_lingui_no2); 

if (retval3 — FALSE) 

{ 

printf(“reture value of close_loop_get_act3 is > or < max, min double\n”); 
return (FALSE); 

> /* End of if retval3 */ 


retval4 - close_loop_get_downup_act4(); 

if (retval4 — FALSE) 

{ 

printf(“reture value of close_loop_get_downup_act4 is > or < max,min doubleW’); 
return (FALSE); 

> /* End of if retval */ 




/* 

* The next for loop calculates the network’s output of the supervise 

* learning 
*/ 
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for (output_lingui_no = 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


f = f2«0.0; 

for (output_term_membership - 0; output_term_membership < 
Output_membership_airaytoutput_lingui_no]; 
output_term_membership++) 

{ 


f +- (Mean_output[output_lingui_no][output_term_membership] * 
Output_term_Deviation[output_lingui_no] [output_term_membership] * 
Act_downupJ4[outputJingui_no][output_term_membership]); 


f2 += (Output_term_Deviation[output_lingui_no][output_term_membership] * 
Act_downup_14[output_lingui_no] [output_term_membership]); 

> /* End of output_term_membership */ 

Control_signal - f / f2; 

if ((Control_signal > MAXFLOAT) II (Control_signal < MINFLOAT)) 

{ 

printf(“Control_signal OUT OF FLOATING POINT RANGED”); 
return (FALSE); 

> 

else 

{ 

printf(“Outputl %1.6f\n”, Control_signal); 

> 

} /* End of output_lingui_no *1 

return (TRUE); 

} 
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#include <stdio.h> 

#include <sx$lib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 
#include“c:\borlandc\file\globalex.h” 


/*********ii^##******#********************************************************»** 

* function close_loop_get_act2() 

* 

* This function calculates the output of the layer two of the network. 

* 

* Input parameter input_lingui_.no, inputx. 

* 

* Return value: int type. 

* 

♦programmer Leung Kam LungMS. Thesis9/28/93 


int 

close_loop_get_act2(unsigned int input_Iingui_no, float inputx) 

{ 

int retval; 


retval - close_loop_activation_12(inputx, input_lingui_no); 
if (retval — FALSE) 

{ 

printf(“reture value of act_function_ptr is > or < max, min double\n”); 
return (FALSE); 

> /* End of if retval */ 

return (TRUE); 

> 

/ 

* *** ************* ******** ********************************* ******* ******* ********** 
*************** 

* This function is called by get_act2 to set up the parameters for calculates the output value 

* for layer of the network. 

* 

♦ Input parameter inputx. 
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*input_lingui_no. 

* 

* Return value: TRUE if successful! calculate the output values. 

* FALSE otherwise. 

* 

* 

* programmer Leung Kam LungMS. Thesis9/28/93 
*/ 

int 

close_loop_activation_12(float inputx, unsigned int input_lingui_no) 

{ 

int membership; 
float temp; 


for (membership - 0; membership < Input_membership_array [input Jingui_.no] ; membership++) 

{ 


temp = Act_12[input_lingui_no][membership] - (float) 
close_loop_12_membership(Mean_input[input_lingui_no][membership), 
Input_term_Deviation[input_lingui_no][membership], 
inputx); 


printf(“a_2 %f temp); 


if (temp >= MAXFLOAT) 

{ 

printf(“close_loop Act_12 value of input_lingui_no %d membership %d is > or < expected 
valueW, 

input_lingui_no, membership); 
return (FALSE); 

> 


> /* End of for membership */ 

printf(‘\n”); 

getchar(); 
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return (TRUE); 


> 


/* End of activation_12 */ 


/*************#****************************************■*■*****■***■**************■*■ 

* The function calculates the activation value for membership function in 

* layer 2 of the NNFLC. 

* input parameters: 

♦mean 
♦deviation 
*u(input value) 

* output: 

♦Return type double 

* 

*programmer:Leung Kam LungMS. Thesis9/28/93 
* 

*/ 


double 

close_loop_I2_inernbership(float mean, float deviation, float u) 

< 


/* 

* printf(“2 deviation %f mean %f float %f ’, deviation, mean, u); 

* getchar(); 

♦/ 

return (exp(-pow(((double) u - (double) mean), 2) / (pow((double) deviation, 2)))); 
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#include <stdio.h> 

#include <stdlib.h> 

#include <conio.h> 

#include “c:\borlandc\file\thesis.h” 
#include“c:\borlandc\file\globalex.h” 


****************************! ****************************************************** 
******* 

* function close_loop_get_act3() 

* 

* This function find the activation value of the layer three by perform the minimun of 

* the output of the input term nodes which are the output of the membership in layer 2. 

* The activation valye of this layer(rule_nodes) are the minimum value from the output 

* of the layer 2. 

* 

* Input parameter input_lingui_nol, 

* input_lingui_no2. 

* 

* 

* Return value: TRUE if successful calculate the output values. 

* FALSE otherwise. 

* 

* 

♦programmerLeung Kam LungMS. Thesis9/28/93 
* 

*/ 

int 

close_loop_get_act3(unsigned int input_lingui_nol, unsigned int input_Iingui_no2) 

{ 

unsigned int rule_node; 

unsigned int lingnol_membership; 
unsigned int lingno2_membership; 

float smallest; 
float smallest2; 

rule_node - 0; 

for (lingnol_membership - 0; lingnoljmembership < MEMBERSHIP_N01; 
lingnol membership++) 

{ 
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smallest - Act_12[input_lingui_nol3[lingnol_membership]; 


for (lingnp2_membership - 0; lingno2_membership < MEMBERSHIP_N02; 
lingno2_membership++) 

{ 


/* printf(“rule_node «%dSt”,rule_node); */ 

smallest2 - Act_12 [input_lingui_no2] [lingno2_membership] ; 

/* printf(“smallest «%f\t”, smallest); *1 
/* printf(“smallest2 ~%f\t”, smallest2); */ 

if (smallest < smallest2) 

{ 

Act_13[rule_node] - smallest; 

/* printf(“Act_13 Act_13[rule_node]); */ 

> 

else 

{ 

Act_13[rule_node] - smallest2; 

/* printf(“Act_13 -%fSn”, Act_13[rule_node]); */ 

} 

rule_node++; 

if (rule_node > RULE_NODE) 

{ 

printf(“ERROR number of rule nodes are greater than the expected valueNn”); 
return (FALSE); 

> 


> /* End of lingno2_membership */ 

> /* End of lingno 1 .membership */ 

return (TRUE); 

> 


I* End of the program */ 



#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 


#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


int 

close_Ioop_get_dovmup_act4(void) 

{ 

unsigned int output Jingui_.no; 
unsigned int output_term_membership; 
unsigned yit rule_index; 
float act4; 

for (output_lingui_no - 0; output_Iingui_no < OUTPUT JLINGUI_NO; 

{ 


for (output_term_membership « 0; output_term_membership < 
Output jrnembership_airay[output_lingui_no]; 
output_term_membership++) 

{ 

act4 - 0.0; 


for (rule_index - 0; rule_index < RULE_NODE; rule_index++) 

{ 


act4 +- Act_13[rule_index] * 

W_14[output Jingui_no] [output_term_membership] [rulejndex] ; 

. > /* End rule_index */ 

/* 

* printf(“ act 4 %8f \n”,act4); getchar(); 

*/ 

if (act4>~ 1) 

{ 

Act_downupJ4[outputJingui_no][output_term_membership] - 1.0; 

> 

else 


t_lingui_no++) 
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{ 


Act_downup_14[output_lingui_no][output_tenn_membership] - act4; 


if (act4 > MAXFLOAT) 

{ 

printf(“activatio value of level 4 is bigger that max floatW’); 
return (FALSE); 

> 

> /* End for output_term_membership */ 


> /* End for output_lingui_no */ 

return (TRUE); 

> 



#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

/ 

**************************** **************************** ************ ************** 
*************** 

* This function find the deviation of each of the membership on layer2. The INPUT_OVERLAP 

* parameter is 2. 

* 

♦Input parameter input_lingui_no 
* 

♦Return value: TRUE if all deviation values are calculated 
♦FALSE otherwise. 

* 

* Programmer:Leung Kam Lung9/30/93MS. Thesis. 

*/ 

int 

find_input_deviation(unsigned int input_lingui_no) 

{ 

float closest_mean; 
float current_mean; 
float temp, 
temp_smallest; 
int membership 1; 
int membership; 

/* printf(“input ling no %d\n”,input_lingui_no); */ 

for (membership = 0; membership < Input_membership_array[input_lingui_no]; membership++) 

{ 


current_mean - Mean_input[input_lingui_no][membership]; 

/* printf(“ current_mean %f\n”,current_mean); */ 

closest_mean - MAXFLOAT; 

for (membershipl - 0; membershipl < Input_membership_array[input_lingui_no]; 
membershipl++) 

{ 

if (membership ! - membership 1 ) 
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{ 


* .* 

temp - Mean_input[input_lingui_no][membershipl]; 

temp_smallest - fabs(temp - current_mean); 

/* printf(“temp smallest %f\n”,temp_smallest); */ 

if (temp_smallest < closest_mean) 

{ 

closest_mean - temp_smallest; 

> 

> 

> /* End of the for membership 1 */ 

if ((closest_mean < MINFLOAT) II (closest_mean > MAXFLOAT)) 

{ 

return (FALSE); 

> 

/* printf(“closest mean %f\n”,closest_mean); */ 

Input_term_Deviation[input_lingui_no] [membership] - (float) (fabs((double) (current_mean - 
closest_mean)) / INPUT_OVERLAP); 

/♦ 

* printf(“Input term D 

* %8f\n”,Input_term_Deviation[input_lingui_no][membership]); 

*/ 

> /* End of the for membership */ 

return (TRUE); 


) /* End of the function find_deviation */ 


/ 

****************************************************** ************** ************** 
*************** 

* This function find the deviation of each of the membership on layer2. The OUTPUT_OVERLAP 

* parameter is 2. 

* 

♦Input parameter input_lingui_no 
* 

♦Return value: TRUE if all deviation values are calculated 
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♦FALSE otherwise. 

* 

' a 

* ProgrammenLeung Kam Lung9/30/93MS. Thesis. 

*/ 

int 

find_output_deviatlon(uiisigned int output_lingui_no) 

{ 

float clo$est_mean; 
float current_mean; 
float temp, 
tempjsmallest; 
int membership 1; 
int membership; 


for (membership - 0; membership < Output_membership_array[output_lingui_no]; membership++) 

< 

current_mean - Mean_output[output_lingui_no][membership]; 

/* printf(“ current_mean %fm”,current_mean); */ 

closest_mean - MAXFLOAT; 

for (membershipl - 0; membershipl < Output_membership_array[output_lingui_no]; 
membership 1 ++) 

{ 


if (membership !- membershipl) 

{ 


temp - Mean_output[output_lingui_no] [membership 1 ]; 

temp_smallest - fabs(temp - current_mean); 

/* printf(“temp smallest %f\n”,temp_smallest); */ 

if (temp_smallest < closest_mean) 

{ 

closest_mean - temp_smallest; 

} 

} 

) /* End of the for membershipl */ 
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if ((closest_mean < MINFLOAT) II (closest_mean > MAXFLOAT)) 

{ •„ 

return (FALSE); 

> 

/* printf(“closest mean %f\n”,closest_mean); */ 

Output_term_Deviation[output_lingui_no][membership] - (float) (fabs((double) (current_mean - 
closest_mean)) / OUTPUTOVERLAP); 

/* 

* printf(“Output term D 

* %8f\n”,Output_term_Deviation[outputJingui_no][membership]); 

*/ 

> /* End of the for membership */ 

return (TRUE); 

> /* End of the function flnd_deviation */ 
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#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include <conio.h> 

#include “c:\borlandcMile\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


/****************************************************************************** 


* This function update the Mean_input of the layer 2 and the Width of the 

* membership function 

* 


* Input parameter Irate, 

* inputx. 

* 

* Output: TRUE if no error otherwise FALSE 

* 


*/ 


int 

error_backpro_layer2(float Irate, int input_index) 

{ 

unsigned int input_lingui_no, 
membership; 

float temp; 

for (input_lingui_no » 0; input_lingui_no < INPUT_LINGUI_NO; input_lingui_no++) 

{ 


. for (membership - 0; membership < Input_membership_array[input_lingui_no]; 
membership-H-) 

{ 


/* clrscr(); */ 

/* 

* printf(‘*\nMEMBERSHIP %d, INPUT_LINGUI_NO %dSn”, membership, 

* input_lingui_no); 

*/ 

Change_E_respect_a - 0.0; 
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rale_cohnection(membership, input_lingui_no); 

/* 

* printf(“lrate %2.5f, index %d, inputx %2.5f\n”, Irate, 

* input_index, Inputx[inputJingui_no][input_index]); 

*/ 

/* 

* printf(“Mean_input %2.5f Diviation %2.5f C_E_respect_a %2.5f 

* act %2.5f\n”, Mean Jnput[inputJingui_no] [membership], 

* Input_term_Deviation[input_lingui_no][membership], 

* Change JE_respect_a, (float) 

* 12_membership(Mean_input[inputJingui_no][membership], 

* Input_term_Deviation[input_lingui_no] [membership] , 

* Inputx[input_Iingui_no][input_index])); 

*/ 

/* Mean_input[input_lingui_no][membership] - */ 

Mean_input_curr[input_lingui_no] [membership] +- temp » 

(Irate * Change_E_respect_a * 

(float) 12_membership(MeanJnput[inputJmgui_no][membership], 
Input_term_Deviation[input_lingui_no][membership], 
Inputx[inputJingui_no][input_index]) * 

(2 * ((Inputx[input_lingui_no][input_index] - Mean_input[input_lingui_no][membership]) / 
(float) pow((double) Input_term_Deviation[input_lingui_no][membership], 2) 

) 

) 

); 


/* printf(“mean curr %2.5f temp); */ 

/* Input_term_Deviation[input_lingui_no] [membership] - */ 

Input_term_Deviation_curr[input_lingui_no] [membership] +- temp - 
(Irate * Change_E_respect_a * 

(float) 12_membership(MeanJnput[input Jingui_.no] [membership], 
Input_term_Deviation[input_linguLno][membership], 
Inputx[input_lingui_no][input_index]) * 

(2 * ((float) pow((double) (Inputx[inputJingui_no][input_index] - 
Mean_input[input_lingui_no][membership]), 2) / 

(float) pow((double) Input_term_Deviation[input_lingui no] [membership], 3) 

) 

) 

); 
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/* printf(“delta deviation %2.5f\n”, temp); */ 

* a 

/* getchar(); */ 


> 

> 

return (TRUE); 

} 


/ 

********************************************************************************** 

****** 

* This function find all connection of rule nodes that are connected to a given 

* in input term membership. 

* 

* Input parameter input_term_membership, 

* input_lingui_.no. 

* 

* Output TRUE if no error otherwise FALSE 
*/ 

int 

rule_connection(unsigned int input_term_membership, unsigned int input_lingui_no) 

{ 

unsigned int N, 
n> 

first_rule_node, 

next_rule_node; 


/* 

* printf(“input_term_membership %d, input_lingui_no %dsn”, 

* input_term_membership, input_lingui_no); 

*/ 

N - (Max_rule_no / (unsigned int) pow((double) MEMBERSHIP_DEMEMSION, 
(double) (INPUT_LINGUI_NO - input Jingui_no))) - 1; 

first_rule_node - (input_term_membership * (unsigned int) pow((double) 
MEMBERSHIP_DEMEMSION, 

(double) (INPUT_LINGUI_NO - input_lingui_no - 1))); 


86 



/* printf (“ mlFIR S T_N ODE %dMi”,first_rule_node); */ 

membership_connection(first_rule_node ) input_lingui_no, input_term_membership); 
find_next_mIe_node(first_rule_node, input_lingui_no, input_term_membership); 


for (n - 1; n <*» N; n++) 

{ 

next_rule_node - first_rule_node + (n * (unsigned int) pow((double) 
MEMBERSHIP.DEMEMSION, 

(double) (INPUTJLINGUI_NO - input_lingui_no))); 

/* printfC^nNEXT_NODE %d \n”, next_rule_node); */ 

membership_connection(next_rule_node, input_lingui_no, input_term_membership); 
find_next_rule_node(next_rule_node > input_lingui_no, input_term_membership); 

> 

/* getcharQ; */ 

return (TRUE); 

> 


/ 

********************************************************************************** 

********* 

* This function find all connection of input term membership that are connected to a given 

* rule node. 

* 

* Input parameter input Jerm_membership, 

* input_node, 

* rule_node. 

* 

* Output TRUE if no error otherwise FALSE 
*/ 
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int 

membership_connection(unsigned int rule_node, unsigned int input_node, 
unsigned int input_term_membership) 


{ 

unsigned int i, 
is_minimum, 
input_lingui_no, 
Rule_no, 
membership; 

float minimum, 
temp_minimum; 


is_minimum « TRUE; 

Rule_no = rule_node; 

if (rule_node >= Max_rule_no li rule_node < 0) 

{ 

printf(“node greater expected value %d’\ rule_node); 
return (FALSE); 

} 

minimum - Act_12[input_node][input_term_membership]; 

/* 

* printf(“input_node %d, input_term_membership %d\n”, input_node, 

* input_term_membership); 

*/ 

for (i - INPUT_LINGUI_NO; i >- 1; i-) 

{ 

input_lingui_no - INPUT_LINGUI_NO - i; 
membership » (rule_node / Rule_connect[i - 1]); 

/* 

* printf(“input_lingui_no %d, membership %d input_lingui_no, 

* membership); 

*/ 

temp_minimum - Act_12[input_lingui_no] [membership]; 
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%■, 4 


/* 

* printf(“minimum %2.5f, temp_minimum %2.5fm”, minimum, 

* temp_minimum); 

*/ 

if (minimum > temp_minimum) 

{ 

is_minimum - FALSE; 

/* printf(“is_minimum is FALSE \n”); */ 
break; 

> 

rule_node %- Rule_connect[i - 1]; 

/* printf(“rule_node now %<Ln”, rule_node); */■ 

} 

if (is_minimum) 

< 

/* 

* printf(“R_no %d, E_respect_a %2.5f Rule_no, 

* Change_E_respect_a); 

*/ 

/* printf(“E_sig_13 %2.5f Error_signal_layer3[Rule_no]); */ 
Change_E_respect_a +- Error_signaI_layer3[Rule_no]; 

/* printf(“Change_E_respect_a %2.5f\n”, Change_E_respect_a); */ 

> 


return (TRUE); 


m 


/ 

**++**+**+*+******************♦***************♦*********************************** 

****** 

* This function find next number of rule nodes that are connected to a given 

* in input term membership, and the current rule node that is calculate. 

* 
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* Input parameter input_term_membership, 

* input_lingui_no, 

*current_rule_node. 

* 

* Output: TRUE if no error otherwise FALSE 
*/ 


int 

find_next_ruIe_node(unsigned int next_node, unsigned int input_lingui_no, 
unsigned int input_term_membership) 


< 

unsigned int loop_node; 

loop_node - (((unsigned int) pow((double) MEMBERSHIP_DEMEMSION, 

(double) (INPUT_LINGUI_NO - input_lingui_no - 1))) - 1); 

/* printf(“ loop node %d \n”,loop_node); */ 

while (loop_node > 0) 

{ 

/* printf(“\nNEXT_NODE %d\n”, next_node + 1); */ 

membership_connection(++next_node, input_lingui_no, input_term_membership); 
loop_node— ; 

> 

return (TRUE); 

> 
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#include <stdio.h> 

#include <stcjlib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


/ 

********************************************************************************** 

************ 

* This function calculate the Error signal for layer 3 

* 

* Input parameter void 

* 

* Output: TRUE if no error otherwise FALSE 

* 

*/ 


int 

error_backpro_layer3(void) 

{ 

unsigned int output_term_membership; 

int output_lingui_no, 
rule_no; 

float temp; 

/* printf(“\nERROR SIGNAL LAYER 3\n”); */ 
for (rule_no = 0; rule_no < RULE_NODE; rule_no++) 

{ 

/* printf(“RULE_NO %d\n”,rule_no); */ 

Error_signal_layer3[rule_no] - 0.0; 

for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 

for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; output_term_membership++) 

{ 

I* printf(“op_t_mship %d, output_term_membership); */ 

/* 
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* printf(“W_14 %2.5f, E_s_14 %2.5f, 

*.iW_14[output_lingui_no][output_tenn_membership][rule_no], 

* Error_signal_layer4[output Jingui_no] [output_term_membershi 

*p]); 

*/ 

if (WJ4[output_lingui_no][output_term_membership][rule_no) > 0.0) 

{ 


Error_signal_layer3[rule_no] += temp = 
Error_signal_layer4[output_lingui_no][output_term_membership]; 

> 

/* 

* printf(“E_s_13 %2.5f, S_E %2.5f\n”,temp, 

* Error_signal_layer3[rule_no]); 

*/ 

> /* End of output_term_membership */ 

} /* End of output_lingui_no */ 

/* getchar(); */ 


} I* End of rule_no */ 

/* getchar(); */ 

return (TRUE); 


} 


/* End of error_backpro_layer3 */ 
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#include <stdio.h> 

#include <st$lib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


/ 

*********************** *********************************************************** 
************ 

* This function calculate the Error signal for layer 4 

* 

* Input parameter index 

* 

* Output: TRUE if no error otherwise FALSE ' 

* 

*/ 

int 

error_backpro_Iayer4(float *Sum_width_act, float *Sum_mean_width_act) 

{ 

unsigned int output_lingui_no, 
i; 

float temp; 

/* 

* Update the Mean_output for each of the output_membership in layer 4 
*/ 


/* printf(“\nNOW IN LAYER 4 V’); */ 

for (output_lingui_no -= 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


/* printf(“Sum_w_act %2.5f Sum_width_acttoutput_lingui_no]); *t 
/* 

* printf(“Sum_m_w_act %2.5fsn”, 

* Sum_mean_width_act[output_lingui_no]); 

*/ 
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for (i - 0; i < Output_membership_array[output_lingui_no]; i++) 

< 

/* printf(“Mean_op %2.5f Mean_output[output_lingui_no][i]); */ 
/* 

* printf(“Out_Devi %2.5f 

* Output_term_Deviation[output_lingui_no] [i]); 

*/ 


Error_signal_layer4[output_lingui_no][i] - temp - (Error_signal_layer5[output_lingui_no] * 

( 

( 

( 

(Mean_output[output_lingui_no][i] * 

Sum_width_act[output_lingui_no] 

)- 

Sum_mean_width_act[output_lingui_no] 

)* 

Output_term_Deviation[output_lingui_no][i] 

)/ 

pow(Sum_width_act[output_lingui_no], 2) 


/* printf(“Error_s_14 %2.5fm”, temp); */ 
> /* End fori*/ 


> /* End for mput_lingui_no */ 

return (TRUE); 

> /* End error_backpro_layrer4 */ 
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#include <stdio.h> 

#include <stdlib.h> 

* A 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


/****************************************************************************** 

* This function calculate the Error of the network and update the Mean_output 

* of the layer 5 and the Width of the membership function 

* 

* Input parameter Irate, 

* index. 

* 

* Output: TRUE if no error otherwise FALSE 

* 

*/ 


int 

error_backpro_layer5(float Irate, unsigned int index) 

{ 

unsigned int output_lingui_no, 
centerjndex, 
devijmdex, 
layer_5, 
i; 

float dev_mean, 
temp, 

sum_dev_mean, 

sum_dev_mean_width; 


/* 

* Update the Mean_output for each of the output_membership in layer 5 
*/ 

dev_mean - 0.0; 

/* printf(“\nNOW IN LAYER 5 Irate %f, index %d\n”,lrate, index); */ 

for (output _lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


Sum_mean_width_act[output_lingui_no] - 0.0; 
Sum_width_act[output_lingui_no] - 0.0; 
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sum_dev_mean - 0.0; 
sum_de;v_mean_width - 0.0; 

for (i - 0; i < Output_membership_array[output_lingui_no]; i++) 
{ 


/* printf(“Act_dup4 %f Act_downup_14[output _lingui_no] [i]) ; */ 

/* 

* printf(“Out_Devi %f 

* Output_terai_Deviation[output_lingui_no][i]); 

*/ 

sum_dev_mean +- temp - (Act_downup_14[output_lingui_no][i] * 
Output_term_Deviation[output_lingui_no][i]); 

/* printf(“sum_dev_meam %2.5f\n”, temp); */ 

> /* End for i */ 

/* printf(“sum_dev_meam %2.5f\n”, sum_dev_mean); */ 

/* getchar(); */ 

Sum_width_act[output_lingui_no] - sum_dev_mean; 

for (center_index - 0; center_index < Output_membership_array[output_lingui_no]; 
center_index++) 

{ 


dev_mean - (Act_downup_14[output_lingui_no][center_index] * 
Output_term_Deviation[output_lingui_no][center_index]); 


/* 

* Calculate the Mean_output of the output layer but not update 

* untill all errors of each layer are calculate. 

*/ 

/* (Mean_output[output_lingui_no][center_index] + */ 
Mean_output_delt[output_lingui_no][center_index] +- temp - 
(Error_array[output_lingui_no] * Irate * 

(dev_mean / sum_dev_mean)); 

/* printf(“dev_mean %2.5f dev_mean); */ 

/* printf(“Mean_output_curr %2.5f temp); */ 
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* printf(“SM_output %2.5f\n”, 

* Mean_output_delt[output_lingui_no][center_index]); 
*/ 


> /* End for center_index */ 

/* getcharO; */ 

for (i - 0; i < Output_membership_array[output_lingui_no]; i++) 

{ 

/* printf(“Mean_op %f Mean_output[output_lingui_no][i]); */ 

/* •’ ' 

* printf(“Out_Devi %f 

* Output_term_Deviation[output_Iingui_no][i]); 

*/ 


/* printf(“Act_dup4 %f Act_downup_14[output_lingui_no][i]); */ 

sum_dev_mean_width +- temp - (Mean_output[output_lingui_no][i] * 
Act_downup_14[output_lingui_no][i] * 
Output_term_Deviation[output_lingui_no][i]); 

/* printf(“sum_d_m_w %2.5f\n temp); */ 

> /* End for i */ 

Sum_mean_width_act[output_lingui_no] - sum_dev_mean_width; 

/* printf(“sum_d_m_w %2.5f\n”, sum_dev_mean_width); */ 

/* getchar(); */ 

/* printf(“Error_a %2.5f\n”, Error_array[output_lingui_no]); */ 

for (devi_index - 0; devijndex < Output_membership_array[output_lingui_no]; 
devi_index++) 

{ 


/* 

* Calculate the Output_term_Deviation of the output layer but 

* not update untill all errors of each layer are calculate. 

*1 
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' A 

/* printf(“devi_i %d devi_index); */ 

/* 

* printf(“Mean_op %2.5f 

* Mean_output[output_linguijno][devi_index]); 

*/ 

/* 

* printf(“Act_dup4 %2.5f ", 

* Act_downup_14[output_lingui_no][devi_index]); 

*/ 

/* Output_term_Deviation[output_lingui_no][devi_index] + */ 

Output_term_Deviation_curr[output_lingui_no][devi_index] +- temp - 

( 

(Error_array[output_lingui_no] * 

Irate 

)• 

( 

( 

( 

(Mean_output[output_lingui_no] [devi_index] * 
sum_dev_mean 
) - sum_dev_mean_width 
)* 

Act_downup_14[output_lingui_no][devi_index] 

)/ 

pow(sum_dev_mean, 2) 

) 

); 

/* 

* printf(“Out_Devi %2.5f, S_out_D %2.5fm”, temp, 

* Output_term_Deviation_curr[output_lingui_no][devi_index]); 

*/ 


> /* End for devi_index */ 


Error_signal_layer5[output_lingui_no] - temp - Error_array[output_lingui_no]; 
/* getcharQ; */ 
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> 


/* End for input_lingui_no */ 


return (TRUE); 

} /* End error_backpro_layrer5 */ 


/ 

********************************************************************************** 

****** 

* This function find all connection of rule nodes that are connected to a given 

* in input term membership. 

* 

* Input parameter input_term_membership, 

* input_lingui_no. 

* 

* Output TRUE if no error otherwise FALSE 
*/ 

int 

rule_connection(unsigned int input_tenn_membership, unsigned int input_lingui_no) 

{ 

unsigned int N, 
n, 

first_rule_node, 

next_rule_node; 


/* 

* printf(“input_term_membership %d, input_lingui_no %d\n”, 

* input_term_membership, input_lingui_no); 

V 

N - (Max_rule_no / (unsigned int) pow((double) MEMBERSHIP_DEMEMSION, 
(double) (INPUT_LINGUI_NO - input_lingui_no))) - 1; 

first_rule_node - (input_term_membership * (unsigned int) pow((double) 
MEMBERSHIP_DEMEMSION, 

(double) (INPUT_LINGUI_NO - input_Iingui_no - 1))); 

printf(‘*mfirst_node %d\n”, first_rule_node); 

membership_connection(first_rule_node, input_lingui_no, input_term_membership); 
find_next_rule_node(first_rule_node, input_lingui_no, input_term_membership); 
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for (n - l;n <- N; n++) 

{ 

next_rule_node - first_rule_node + (n * (unsigned int) pow((double) 

MEMBERSHIP JDEMEMSION, 

(double) (INPUT_LINGUI_NO - input_lingui_no))); 

printf(“Snnext_node %d\n”, next_rule_node); 

membership_connection(next_rule_node, input_lingui_no, input_term_membership); 
find_next_rule_node(next_rule_node, input_Iingui_no, input_term_membership); 

> 

getchar(); 
return (TRUE); 

> 


/ 

********************************************************************************** 

3#eafe3fc3|ea#cafe3|ea|e3fe 

* This function find all connection of input term membership that are connected to a given 

* rule node. 

* 

* Input parameter input_term_membership, 

* input_node, 

* rule_node. 

* 

* Output: TRUE if no error otherwise FALSE 
*/ 


int 

membership_connection(unsigned int rule_node, unsigned int input_node, 
unsigned int input_term_membership) 


{ 

unsigned int i, 
is_minimum, 
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input_lingui_no, 

Rule_nQ, 

membership; 

float minimum, 
temp_minimum; 

is_minimum - TRUE; 

Rule_no = rule_node; 

if (rule_node >-= Max_rule_no II rule_node < 0) 

{ 

printf(“node greater expected value %d”, rule_node); 
return (FALSE); 

> 

minimum - Act_12iinput_node][input_term_membership]; 

/* 

* printf(“input_node %d, input_term_membership %dSn”, input_node, 

* input_term_membership); 

*/ 

for (i - INPUT_LINGUI_NO; i >- 1; i-) 

{ 

input_lingui_no - INPUT_LINGUI_N O - i; 
membership - (rule_node / Rule_connect[i - 1]); 

printf(“input_lingui_no %d, membership %<Fn”, input_lingui_no, membership); 

temp_minimum - Act_12[input_lingui_no][membership]; 

printf(“minimum %2.5f, temp_minimum %2.5fSn”, minimum, temp_minimum); 

if (minimum > tempjminimum) 

{ 

is_minimum ** FALSE; 
printf(“is minimum is FALSE Nn”); 
break; 

> 
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rule_node %= Rule_connect[i - 1]; 

/* prii}tf(“rale_node now %d\n”, rule_node); */ 

> 

if (is_minimum) 

{ 

printf(“R_no %d, E_respect_a %2.5f Rule_no, Change_E_respect_a); 
printf(“Error_signal_layer3 %2.5f Eiror_signal_layer3[Rule_no]); 
Change_E_respect_a += Error_signal_layer3[Rule_no]; 
printf(“E_respect_a %2.5f\n”, Change_E_respect_a); 

> 

return (TRUE); 


> 

/ 

********************************************************************************** 

****** 

* This function find next number of rule nodes that are connected to a given 

* in input term membership, and the current rule node that is calculate. 

* 

* Input parameter input_term_membership, 

* input_lingui_no, 

*cuirent_rule_node. 

* 

* Output: TRUE if no error otherwise FALSE 
*/ 


int 

find_next_ru!e_node(unsigned int next_node, unsigned int input_lingui_no, 
unsigned int input_term_membership) 


{ 

unsigned int loop_node; 

loop_node - (((unsigned int) pow((double) MEMBERSHIP J3EMEMSION, 
(double) (INPUTJJNGUI_NO - input_lingui_no - 1))) - 1); 

/* printf(“ loop node %d\n”,loop_node); */ 
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while (loop_node > 0) 

( 

printf(‘*\nnext node %d\n”, next_node + 1); 

membership_connection(++next_node, input_lingui_.no, input_term_membership); 
loop_node— ; 

> 

return (TRUE); 

> 
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#include <stdio.h> 

#include <stfilib.h> 

#include <conio.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


/********* ******** ************************************ ********************** ******* 

* function updata_input_weight(). 

★ 

* This function updata the weight of the membership function which has the closest 

* “mean” for a input data(inputx) of an input linguistic node identify as 

* input_ling_index. The find_mini_ptr return the integer indicating the index 

* of the membership function. 

* This function call the call the function fi_mini_index() by referrencing the 

* function address thus the find_mini_ptr is pointed to function fi_mini_index 

* 

* Input parameter 
*inputx, 

*alpha, 

*input_ling_index, 

* (* find_mini_ptr)() 

* 

* Output: NONE 

* 

* Programmer 

* Leung Kam Lung 9/14/93MS. Thesis. 

* 

*/ 

void 

update_input_12_w(float inputx, float alpha, unsigned int input_ling_node) 

{ 


unsigned int mini_index; 

unsigned int (*find_mini_ptr) (float, unsigned int); 

find_mini_ptr - fi_12_mini_index; 

mini_index - (*fmd_mini_ptr) (inputx, input_ling_node); 

/* 

* printf(“input node %d mini_index %d before %f “,input_ling_node, 

* mini_index, Mean_input[input_ling_node][mini_index]); 

*/ 
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Mean_input[input_ling_node][mini_index] = (Mean Jnput[inputJing_node] [mini Jndex] + (alpha * 
. a (inputx - Mean_input[input_ling_node][mini_index]))); 

/* printf(“after %Ni”,Mean_input[input_ling_node][mini_index]); */ 


> 

/************#**************************************************************** 

* function fi_mini_index(); 

* 

* This function find the shortest distance between the inputx and the “mean” 

* of each of the membership function in a set of membership function which 

* are belong to a particular input linguistic node. 

* 

* Input parameter 
*inputx, 

*input_ling_node. 

* 

* Return value: mini_index 

* 

* programmerLeung kam Lung9/16/93MS. Thesis. 

* 

*/ 


unsigned int 

fi_12_mini_index(float inputx, unsigned int input_ling_node) 

{ 

double minimum, 
temp, 

distance[MAX_MEMBERSHIP]; 

unsigned int k, 
mini_index, 
membership; 

k-0; 
temp - 0.0; 
mimjndex - 0; 
membership - 0; 
minimum - 0.0; 

/* printf(“input_ling_node %dSn”,input_ling_node); */ 
k - Input_membership_array[input_ling_node]; /* k contain the number of 

* membership function for 

* the input_ling_node */ 
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for (membership - 0; membership < k; membership++) 

{ 


' / 

* calculates the distance between the inputx and the “meam” of each 

* membership function 
*/ 

distancefmembership] - pow((double) (inputx - Mean_input[input_ling_node][membership]), 2); 

> 


/* 

* Find the shortest distance between inputx and “meam” of all membership 

* functions for a given input linguistic node. 

*/ 

minimum « distance[0]; 

/* printf(“minimum %fn”, minimum); */ 

for (membership = 1; membership < k; membership++) 

{ 

temp » distance[membership]; 

/* printf(“temp minimum %fn”,temp); */ 

if (temp < minimum) 

{ 

minimum - temp; 

/* printf(“minimum %f\n”,minimum); */ 
mimjndex - membership; 

/* printf(“ minimum index %d^n”,mini_index); */ 

> 

} 

return (mimjndex); 

> 


/************************** ******* *********************** ♦♦*j)t*J*C3*t*%)tt5ft***Jtc**3jC*j(C** **** 

* function updata_output_w(). 

* 

* This function updata the weight of the membership function which has the closest 
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* “mean” for a output data(inputy) of an output linguistic node identify as 

* output_ling_index. The fmd_output_mini_ptr return the integer indicating the index 

* of the membership function. 

* This function call the call the function fi_out_mini_index() by referrencing the 

* function address thus the find_output_mini _ptr is pointed to function fi_out_mini_index 

* 

* Input parameter 
♦inputy, 

♦alpha, 

*output_ling_index, 

* (* find_output_mini_ptr)() 

* 

* Output: NONE 

* 

* Programmer 

* Leung Kam Lung 9/14/93MS. Thesis. 

* 

*/ 

void 

update_output_w(float inputy, float alpha, unsigned int output_ling_node) 

{ 

unsigned int mini_index; 

unsigned int (*fmd_output_mini_ptr) (float, unsigned int); 

find_output_mini_ptr ** fi_output_mini_index; 

mini_index - (*find_output_mini_ptr) (inputy, output_ling_node); 

/* 

* printf(“output node %d mini_index %d before %f “,output_ling_node, 

* mini_index, Mean_output[output_ling_node][mini_index]); 

*/ 

Mean_output[output_ling_node][mini_Lndex] - (Mean_output[output_ling_node][mini_index] + 
(alpha * (inputy - Mean_output[output_ling_node][mini_index]))); 


/* 

* printf(“ Mean_output - 

* %f\n”,Mean_output[output_ling_node] [mini_index]); 

*/ 

/* printf(“after %f\n”,Mean_output[output_ling_node][mini_index]); */ 


> 
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/*************** *************** *********************************************** 

* function fi_output_mini_index(); 

* 

* This function find the shortest distance between the inputy and the “mean” 

* of each of the membership function in a set of membership function which 

* are belong to a particular output linguistic node. 

* 

* Input parameter 
♦inputy, 

*output_ling_node. 

* 

* Return value: mini_index 

* 

* programmerLeung kam Lung9/16/93MS. Thesis. 

* 

*/ 

unsigned int 

fi_output_mini_index(float inputy, unsigned int output_iing_node) 

{ 

double minimum, 
temp, 

distancetMAX.MEMBERSHIP]; 

unsigned int k, 
minLindex, 
membership; 

k “ 0; 
temp - 0; 
mini_index - 0; 
membership - 0; 
minimum = 0.0; 



k - Output_membership_array[output_ling_node]; /* k contain the number 

* of membership function 

* for the 

* input_ling_node */ 

for (membership - 0; membership < k; membership++) 

{ 

/* 

* calculates the distance between the inputx and the “meam” of each 
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* membership function 
*/ 

distance[membership] - pow((double) (inputy - Mean_output[output_ling_node][membership]), 

2 ); 

> 

7 * 

* Find the shortest distance between inputx and “meam” of all membership 

* functions for a given input linguistic node. 

*/ 

minimum - distance[0]; 

for (membership - 1; membership < k; membership++) 

{ 

temp - distance[membership]; 
if (temp < minimum) 

{ 

minimum - temp; 
minijndex - membership; 

> 

> 

return (minijndex); 

} 
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#include <stdio.h> 

#include <stdlib.h> 

#include <conio.h> 

#include <math.h> 

#include “c:\borlandc\file\thesisir 
#include “c:\borlandc\file\globalex-h” 

/ 

********************************************************************************** 

********* 

* This function get the mean value of all the membership function for all input linguistic 

* nodes and all the membership function for all output linguistic nodes. 

* 

* Input parameter Global paramters; 

*Max_range, Min_range, 

* Membership_airay, Mean_inpuL 

* 

* Output 

♦TRUE :If all mean values were calculted. 

♦False: If any mean values were not calculted. 

* 

* programmenLeung Kam LungMS. Thesis9/28/93. 

* 

* 

♦All input and output variables are array pointer. 

* 

*/ 

unsigned int 
get_means(void) 

{ 

unsigned int input_no, 
output_no, 
membership_index; 

clrscr(); 

for (input_no - 0; input_no < INPUT_LINGUI_NO; input_no++) 

{ 


if (find_mean(Max_range[input_no], Min_range[input_no], 
Input_membership array[input no], input_no)) 

{ 


for (membership_index - 0; membership_index < Input_membership_array[input_no]; 
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membership_index++) 

{ 

* A 

printf(“Initail Mean_input%d - %f\n”, membership_index, 
Mean_input[input_no][membership_index]); 

> /* End of for membership Jndex */ 

> 

else 

{ 

printf(“fail to find input meanW); 
return (FALSE); 

> 

getchar(); 

> /* End of for input_.no */ 

for (output_.no - 0; output_.no < OUTPUT_LINGUI_NO; output_no++) 

{ 


if (find_output_mean(Output_max_range[output_no], Output_min_range[output_no], 
Output_membership_array[output_no],output_no)) 


{ 


for (membershipjndex - 0; membership_index < Output_membership_array[output_no); 
rnembership_mdex++) 

{ 


printf(“Initiall Mean_output%d - %f\n”, membership_index, 
Mean_output[output_no][membership_index]); 

} /* End for membership_index */ 

> 

else 

{ 


printf(“fail to find output meanNn”); 
return (FALSE); 


> 

getchar(); 

) /* End of for output_.no */ 

return (TRUE); 


111 



> 


/************^***** ******************************************************** ******* 

* This function calculate the mean value for input values of range from 

* MAX_RANGE to MIN_RANGE. The mean value is start from left to right for 

* a set of input term nodes which are correspond to a input linguistic 

♦node. 

* 

* input parameters 

* MAX_RANGE 

* MIM_RANGE 

* NUM_MEMBERSHIP 

* MEAN 

* 

* Output: TRUE if no arithmatic error 

* FALSe otherwise. 

* 

* programmenLeung Kam LungMS. Thesis9/28/93. 

*/ 


unsigned int 

find_mean(float max range, float min_range, unsigned int num_membership, unsigned int 
input_no) 

{ 


float segment; 
unsigned int i; 


/* 

* The value of segment descipt equal distance between each mean. 

*/ 

printf(“max_range %f , min_range %f, inputjnoW’, max_range, min_range, input_no); 
segment - (max_range - min_range) / (num_membership +1); 

if ((segment > MAXFLOAT) II (segment < MINFLOAT)) 

{ 

return (FALSE); 

> 

getchar(); 

/* 
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* The value of i start from NUM__MEMBERSHIP and decrease to 1 

* 

*/ 

Mean_input[input_no][0] - min_range + segment; 
for (i - 1; i < num_membership; i++) 

< 

Mean_input[input_no][i] - Mean_input[input_no][i - 1] + segment; 

} 

return (TRUE); 

> 


/******************** ****************^******************************************** 

* This function calculate the mean value for output values of range from 

* MAX_RANGE to MIN_RANGE. The mean value is start from left to right for 

* a set of output term nodes which are correspond to a output linguistic 

* node. 

* 

* input parameters 

* MAX.RANGE 

* MIM_RANGE 

* NUM_MEMBERSHIP 

* MEAN 

* 

* Output TRUE if no arithmatic error 

* FALSe otherwise. 

* 

* programmenLeung Kam LungMS. Thesis9/28/93. 

*/ ~ 

unsigned int 

find_output_mean(float max_range, float min_range, unsigned int num_membership, unsigned int 
output_.no) 

{ 


float segment; 
unsigned int i; 


/* 

* The value of segment descipt equal distance between each mean. 
*/ 


segment - (max_range - min_range) / (num_membership + 1); 



/* printf(“segment %f\n”, segment); */ 

* A 

if ((segment > MAXFLOAT) II (segment < MINFLOAT)) 

{ 

return (FALSE); 

> 

/* 

* The value of i start from NUM_MEMBERSHIP and decrease to 1 

* 

*/ 

Mean_output[output_no][0] - minjrange + segment; 

for (i - 1; i < num_membership; i++) 

{ 

Mean_output[output_no][i] = Mean_output[output_no][i - 1] + segment; 

> 

return (TRUE); 

> 
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sassaaaagsaaasa 

wmwmmwmwwwwwwww 


#include “c:\borlandc\file\thesis.h” 


FILE *Inputx_fptr, 

*Output_fptr, 

*Input_term_Deviation_fptr, 

*Final_input_term_deviation_fptr, 

*Output_term_Deviation_fptr; 

*Final_output_term_deviation_fptr, 

*Mean_output_fptr, 

*Final_mean_output_fptr, 

*Mean_input_fptr, 

*Final_mean_input_fptr; 

*W_14_fptr, 

*Error_fptr, 

*Fback_fptr, 

*Consigna_fptr, 

*Slope_fptr, 

*Max_data_fptr, 

int Rule_connect[INPUT_LINGUI_NO]; 

/* 

* These values are for filtering the input and output data. 

*/ 

float Max_inputl, 

Max_input2, 

Max_outputl; 
float Min_inputl, 

Min_input2, 

Min_outputl; 
float Input_slopel, 

Input_slope2, 

Control_signal, 

Output_slopel; 

float Max_range[INPUT_LINGUI_NO]; /* input range of linguistic nodes */ 
float Min_range[INPUT_LINGUI_NO]; /* input range of linguistic nodes */ 
float Output_max_range[OUTPUT_LINGUI_NO]; /* output range of linguistic nodes */ 
float Output_min_range[OUTPUT_LINGUI_NO]; /* output range of linguistic nodes */ 
float Input_term_Deviation[INPUT_LINGUI_NO] [M AX_MEMBERSHIP] ; /* input term deviation */ 
float Input_term_Deviation_curr[INPUT_LINGlJI_NO][MAX_MEMBERSHIP]; /* input term devia- 
tion */ 

float Output_term_Deviation[OUTPUT_LINGUI_NO][MAX_MEMBERSHIP]; /* output term devia- 
tion */ 
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float Output_term_Deviation_curr[OUTPUT_LINGUI_NO] [MAX_MEMBERSHIP] ; /* output term 
deviation */, 

float Inputx[INPUT_LINGUI_NO] [MAX_INPIJT_INDEX] ; /* trainning data */ 
float Output[OUTPUT_LINGUI_NO] [MAX_INPUT_INDEX] ; /* trainning data */ 
float Sum_width_act{OUTPUT_LINGUI_NO]; /* sum of width and activation of 
* layer */ 

float Sum_mean_width_act[OUTPUT_LINGUI_NO]; /* sum of width and activation 

* and mean of layer */ 

float Segment; /* number of segment of the range MAXIMIN 

* ranger */ 


/* 

* mean value of each membership of output term nodes 
*/ 


float Mean_output[OUTPUT_LINGUI_NO][MAX_MEMBERSHIP]; 
float Mean_output_delt[OUTPUT_LINGUI_NO][MAX_MEMBERSHIP]; 


/* 

* mean value of each membership of input term nodes 
*/ 


float Mean_input[INPUT_LINGUI_NO][MAX_MEMBERSHIP]; 
float Mean_input_curr[INPUT_LINGUI_NO][MAX_MEMBERS.HIP]; 

float Close_loop_input[INPUT_LINGUI_NO] ; 


/* input number use for close loop simulation */ 
float Cal_enror, /* calculate output error */ 

/* float Output; */ 

float ReffMAXJTIME]; 


/* 

* weight of input term nodes for all input linguistic nodes 
*/ 

float WJ2[MAX_MEMBERSHIP][INPUT_LINGUI_NO]; 

/* 

* activation value of layer 2 from input linguistic node 1 and node 2 
*/ 
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float Act_12[INPUT_LINGUI_NO] [MAX_MEMBERSHIP]; 

' *» 

/* 

* weight in layer 3 
*/ 

float W_13 [RULE_NODE] [ALL_L2_MEMBERSHIP]; 

/* 

* activation value of each rule nodes in layer 3 from each input term nodes of layer 2 
*/ 

float Act_13 [RULE_NODE] ; 

/* 

* error signal of layer 3, error for each of the rule node. 

*/ 

float Error_signal_layer3[RULE_NODE]; 

float W_14[OUTPUT_LINGUI_NO][OUTPUT_TERM_MEMBERSHIP][RULE_NODE]; /* weight 
in layer 4 */ 

/* 

* activation value of each output term nodes in layer 4 from each rule nodes of layer 3 
*/ 

float Act_14[OUTPUT_LINGUI_NO][OUTPUT_TERM_MEMBERSHIP]; 
float Act_downup_14[OUTPUT_LINGUI_NO][OUTPUT_TERM_MEMBERSHIP]; 


float Output_up_down; 
float Change_E_respect_a; 

float Output_down_up[OUTPUT_LINGUI_NO][MAX_INPUT_INDEX]; 

int Max_rule_no; 
int Before_loop; 
int Max_indexs; 

unsigned int INPUT_NODE; 
unsigned int OUTPUT_NODE; 

/* 

* contain the membership value of each of the input linguistic input node 
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*/ 

* i» 

unsigned int Input_membership_array[INPUT_LINGUI_NO]; 

unsigned int Output_membership_array[OUTPUT_LINGUI_NO]; 
float Error_array[OUTPUT_LINGUI_NO]|>IAX_INPUT_INDEX]; 
float Error_signal_layer5 [OUTPUT_LINGUI_NO] ; 

float Error_signaljayer4[OUTPUT_LlNGUI_NO][OUTPUT_TERM_MEMBERSHIP]; 
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#include <stdio.h> 

#include <s.tdlib.h> 

#include <conio.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

unsigned 

initials_aIl(void) 

< 

unsigned int output_lingui_no, 
center_index, 
rule_no, 
i, 

devijndex; 

init_maxmin_range(); 

get_means(); 

init_weight2(); 

init_weight3(); 

init_weight4(); 

for (output lingui_no - 0; output Jingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 

for (centerjndex - 0; center_index < Output_membership_array[output_lingui_no]; 
center_index++) 

{ 

Mean_output_delt[outputJingui_no] [center Jndex] - 0.0; 

> 

for (devi_index = 0; devijndex < Output_membership_array[output_lingui_no]; 
devi_index++) 

{ 

Output_term_Deviation_curr[output_lingui_no] [devijndex] - 0.0; 

> 

> 

Rule_connect[0] - 1; 

Rule_connect[l] - lnput_membership_array[0]; 

for (i - 2; i < INPUT_LINGUI_NO; i++) 

{ 

Rule_connect[i] - (Input_membership_array[i - 1] * Rule_connect[i - 1]); 

> 
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if (INPUT_LINGUI_NO > 1) 

{ 

Max_rule_no - Rule_connect[INPUT_LINGlJI_NO - 1] * 
Input_membership_array[INPUT_LINGUI_NO - 1]; 

> 

printf(“ Max rule no %d\n”, Max_rule_no); 

getchar(); 

return (TRUE); 

} 
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#include “c:\borIandc\file\thesis.h” 

float Max_range [INPUT_LINGUI_NO] ; /* range of linguistic nodes */ 
float Min_range[INPUT_LINGUI_NO]; /* range of linguistic nodes */ 
float Segment; /* number of segment of the range MAXIMIN 

* ranger */ 


I* 

* mean value of each membership of output term nodes 
*/ 


float Meam_outputtOUTPUT_TERM_MEMBERSHIP]; 

/* 

* mean value of each membership of input term nodes 
*/ 

float Mean_input[INPUT_LINGUI_NO] [MAX_MEMBERSHIP] ; 

/* 

* weight of input term nodes for all input linguistic nodes 
*/ 

float W_12[MAX_MEMBERSHIP][INPUT_LINGUI_N0]; 

/* 

* activation value of layer 2 from input linguistic node 1 and node 2 
*/ 

float Act_12[INPUT_LINGUI_NO][MAX_MEMBERSHIP]; 

/* 

* weight in layer 3 
*/ 

float W_13 [RULE_NODE] [ALL_L2_MEMBERSHIP]; 

/* 

* activation value of each rule nodes in layer 3 from each input term nodes of layer 2 
*/ 

float Act_13[RULE_NODE]; 

float WJ4[OUTPUT_TERM_NO][RULE_NODE]; /* weight in layer 4 */ 

/* 
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* activation value of each output term nodes in layer 4 from each rule nodes of layer 3 
*/ 

float Act_14[OUTPUT_TERM_NO]; 
float Output_up_down; 

float W_15_output_up_down[OUTPUT_LINGUI_NO]; 

unsigned int INPUT_NODE; 
unsigned int OUTPUT_NODE; 

/* 

* contain the membership value of each of the input linguistic input node 
*/ 

unsigned int Membership_array[INPUT_LINGUI_NO]; 


void 

initgloble(void) 

{ 

Membership_array[0] - MEMBERSHIP_N01; 

Membership_anay[l] - MEMBERSHIP_N02; 

Max_range[0] - 10.0; 

Min_range[0] - -10.0; 

Max_range[l] - 10.0; 

Min_range[l] - -10.0; 

NUM.MEMBERSHIP - 0; 

> 

void 

init_weight(void) 

{ 

for (j - 0; j < Membership_array[input_lingui_no]; j++) 

{ 

for (input_lingui_.no - 0; input_lingui_no < INPUT_LINGUI_NO; input_lingui_no++) 

{ 

WJ2[j] [input_lingui_no] - random(); 

> 

> 

> 
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#include <time.h> 

#include “c: N Jborlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


void 

init_maxmin_range(void) 

{ 

unsigned int input_lingui_no, 
input_index, 
output_lingui_no, 
output_index; 

float temp, 
temp_min, 
temp_max; 

Input_membership_array [0] - MEMBERSHIP_NO 1 ; 
Input_membership_airay[l] - MEMBERSHIP_N02; 
Output_membership_array[0] - OUTPUT_TERM_MEMBERSHIP; 
temp - 0; 

for (input_lingui_no - 0; input_lingui_no < INPUT_LINGUI_NO; 
input_lingui_no++) 

< 

temp_min - lnputx[input_lingui_no][0]; 
temp_max - lnputx[input_lingui_no][0]; 

for (input_index - 1; input_index < MAX_INPUT_INDEX; input_index++) 

{ 


temp - Inputx[input_lingui_no][input_index]; 

/* printf(“ temp %f\n”,temp); */ 

if (temp < temp_min) 

{ 

temp_min - temp; 

/* 

* printf(“ temp_min %fm”,temp_min); getcharO; 
*/ 

} 

else 

{ 
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if (temp > temp_max) 
temp_max - temp; 

/* 

* printf(“ temp_max %f\n”,temp_max); getcharO; 
*/ 



> 


/* 

* printf(“max %f, min %f, input_no %cf\n”, temp_max, temp_min, 

* input_lingui_no); 

*/ 

Max_range[input_Iingui_no] - temp_max; 
Min_range[input_iingui_no] - temp_min; 

/* 

* printf(“ max ramge of input_lingui_no %d is %f\n”,inputJingui_nQ, 

* temp_max); printf(“ min ramge of input_lingui_no %d is 

* %f\n”,input_lingui_no, temp_min); getchar(); 

*/ 

temp_min - 0; 
temp_max - 0; 


> /* End of inpot_lingui_.no */ 

temp - 0; 

for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGlJI_NO; 
output_Iingui_no++) 

{ 


temp_min - Output[output_lingui_no][0]; 
temp_max - Output[output_lingui_no][0]; 

for (outputjndex - 1; outputjndex < MAX INPUTJNDEX; output_index++) 

{ 


temp - Output[outputJingui_no] [outputjndex]; 

if (temp < temp_min) 

{ 

tempjmin - temp; 

> 
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else 

< ,, 

if (temp > temp_max) 

{ 

temp_max - temp; 

} 

> 

} /* End of outputjndex */ 


Output_max_range[output_lingui_no] - temp_max; 
Output_min_range [output_lingui_.no] - temp_min; 

/* 

* printf(“ max ramge of output_lingui_no %d is 

* %f\n’\output_lingui_no, temp_max); printf(“ min ramge of 

* output_lingui_no %d is %f\n”, output_lingui_no, temp_min); 

* getchar(); 

*/ 

tempjmin - 0; 
temp_max - 0; 


} /* End of output_lingui_.no */ 


> 

unsigned int 
init_weight2(void) 
{ 


unsigned int input_lingui_no, 
membership_index; 

randomizeO; 

for (input_lingui_no - 0; input Jingui_no < INPUT_LINGUI_NO; 
input_lingui_no++) 

{ 

for (membership_index - 0; membership Jndex < Input_membership_array[inputJingui_no]; 
membership_index++) 

{ 

WJ2[input_lingui_no][membership_index] - (random( 10001) 1 10000.0); 

/* 

* printf(“input_lingui_no %d membership_index%d, weight - 
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* %3.5f\n”, input_lingui_no, membership_index, 

* WJ2[input_lingui_no][membershipJndex]); 

*/ 


> 

> 

return (TRUE); 

> 


unsigned int 
init_weight3(void) 
{ 


unsigned int rule_node_index, 
all_12_membership_index; 

randomize(); 

for (rule_node_index - 0; rule_node_index < RULE_NODE; 
rule_node_index++) 

{ 

for (all_12_membership_index - 0; all_12_membership_index < ALL_L2_MEMBERSHIP; 
all_12_membership_index++) 

{ 

W_13 [rule_node_index] [all _12_membership_index] - (random( 10001) / 10000.0); 

/* 

* printf(“rule_node %d all_12_membership%d, weight - %3.5fn”, 

* rule_node_index, all_12_membership_index, 

* W_13[rule_node_index][all_12_raembership_index]); 

*/ 


} 

> 

return (TRUE); 

> 


unsigned int 
init_weight4(void) 
{ 


unsigned int output_term_node_index, 
rule_node_index, 
output_lingui_node; 
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for (output_lingui_node - 0; output_lingui_node < OUTPUT_LINGUI_NO; output_lingui_node++) 

< 

for (output_term_node_index - 0; output_term_ncxie_index < OUTPUT_TERM_MEMBERSHIP; 
output_term_node_index++) 

{ 

for (rule_node_index - 0; rule_node_index < RULE_NODE; 
rule_node_index++) 

{ 

W_14[output_lingui_node] [output_term_node_index] [rule_node_index] ** 1.0; 

/* 

* printf(“output_term_node_index %d, rule_node %d, weight - 

* %3.5f\n”, output_term_node_index, rule_node_index, 

* W_14[output_term_node_index] [rule_node_index]); 

* 

*/ 

} 

> 

> 

return (TRUE); 

> 
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finclude <stdio.h> 

#include <$tdlib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex.h” 


int 

get_act2(void) 

{ 

int retval; 

unsigned int input_lingui_no; 
unsigned int input_index; 


int (*act_function_ptr) (float, unsigned int); 
act_function_ptr - activation_12; 

for (input_lingui_no - 0; input_lingui no < INPUT_LINGUI_NO; input_lingui_no++) 

{ 


for (input_index - 0; input_index < MAXJNPUTJNDEX; input_index++) 

{ 


retval - (*act_function_ptr) (Inputx[input_lingui_no][input_index], input_lingui_no); 


} 


if (retval — FALSE) 

{ 


printf(“reture value of act_function_ptr is > or < max, min doubleNn”); 
return (FALSE); 

} /* End of if retval */ 

> /* End of for input_index */ 

> /* End of for input_lingui_no */ 

return (TRUE); 


int 

activation_12(float inputx, unsigned int input_lingui_no) 

{ 

int membership; 
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float temp; 


double (*activation_12_ptr) (float, float, float); 
activation_12_ptr - 12„membership; 

for (membership - 0; membership < Input_membership_array[input_lingui_no]; membership++) 

{ 


temp - Act_12[mput_lmgui_no] [membership] - (float) (*activation_12_ptr) 
. (Mean_input[input_lingui_no] [membership], 

Input_term_Deviation[input_lingui_no][membership], 

inputx); 


if ((temp > MAXFLOAT) II (temp < MINFLOAT)) 

{ 


printf(“Act_12 value of input_lingui_no %d membership %d is > or < expected valueNn”, 
input_lingui_no, membership); 

return (FALSE); 

> 


> /* End of for membership */ 

return (TRUE); 


> /* End of activation_12 */ 


j * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * $ * * * * * * * * # # ifc 

* The function calculates die activation value for membership function in 

* layer 2 of the NNFLC. 

* input parameters: 

♦mean 
♦deviation 
*u(input value) 

* output: 

♦return type double 

* 
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*/ 


double 

I2_membership(fioat mean, float deviation, float u) 

{ 


return (exp(-pow(((double) u - (double) mean), 2) / 
(2 * pow((double) deviation, 2))) / 
(2.50662825 * (double) deviation)); 

> 
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#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include “c:\Jborlandc\file\lhesis.h” 

#include “c:\borlandc\file\globalex.h” 

/******************************************************************************* 

* function get_act2() 

* 

* This function calculates the output of the layer two of the network. This 

* function uses a function pointers points to function activation_12. 

* 

* Input parameter none. 

* 

* Return value: int type. 

* 

*programmer:Leung Kam LungMS. Thesis9/28/93 
* 

*/ 


int 

get_act2(unsigned int inputlinguino, unsigned int input_index) 

{ 

int retval; 

int (*act_function_ptr) (float, unsigned int); 
act_function_ptr - activation_12; 

retval - (*act_function_ptr) (Inputx[input_lingui_no][input_index], input_lingui_no); 
if (retval — FALSE) 

{ 

printf(“ieture value of act_function_ptr is > or < max, min double\n”); 
return (FALSE); 

} /* End of if retval */ 

return (TRUE); 

> 


*********************************************************************** 4 *********** 
*************** 

* This function is called by get_act2 to set up the parameters for calculates the output value 
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* for layer of the network. 

* 

' a 

* Input parameter inputx, 

*input_lingui_no. 

* 

* Return value: TRUE if successful calculate the output values. 

* FALSE otherwise. 

* 

* 

♦programmenLeung Kam LungMS. Thesis9/28/93 
*/ 

int 

activation_12(float inputx, unsigned int input_Iingui_no) 

< 

int membership; 
float temp; 

double ( * acti vation_12_ptr) (float, float, float); 
activation_12_ptr - 12_membership; 

for (membership - 0; membership < Input_membership_array[input_lingui_no]; membership++) 

{ 


temp = Act_12[input_lingui_no][membership] - (float) (*activation_12_ptr) 
(Mean_input[input_lingui_no] [membership] , 
Input_term_Deviation[input_lingui_no][membership], 
inputx); 


/* 

* printf(“act layer2 %f\n”,temp); getchar(); 
*/ 


if (temp >- MAXFLOAT) 

{ 


printf(“Act_12 value of input _lingui_.no %d membership %d is > or < expected valueV’, 
input_lingui_no, membership); 

return (FALSE); 

> 
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> 


/* End of for membership ♦/ 


return (TRUE); 


> 


/♦ End of activation_12 */ 


/*********************************************** ********************* ********** 

* The function calculates the activation value for membership function in 

* layer 2 of the NNFLC. 

* input parameters: 

♦mean 
♦deviation 
*u(input value) 

* output: 

♦Return type double 

* 

♦programmenLeung Kam LungMS. Thesis9/28/93 
* 

*/ 


double 

12_membership(float mean, float deviation, float u) 


/♦ 

* printf(“2 deviation %f mean %f float %f ’, deviation, mean, u); 

* getchar(); 

*/ 

return (exp(-pow(((double) u - (double) mean), 2) / (pow((double) deviation, 2)))); 

> 



#include <stdio.h> 
#include <stdlib.h> 
#include <conio.h> 


#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

/ 

************************************** ********************* ************ *********** 
******* 

* function get_act3() 

* 

* This function find the activation value of the layer three by perform the minimun of 

* the output of the input term nodes which are the output of the membership in layer 2. 

* The activation value of this layer(rule_nodes) are the minimum value from the output 

* of the layer 2. 

* 

* Input parameter input_lingui_noI, 

* input_lingui_no2. 

* 

* 

* Return value: TRUE if successful! calculate the output values. 

* FALSE otherwise. 

* 

* 

♦programmer: Leung Kam LungMS. Thesis9/28/93 
* 

*/ 

int 

get_act3(unsigned int input_lingui_nol, unsigned int input_Iingui_no2) 

{ 

unsigned int rule_node; 

unsigned int lingnol_membership; 
unsigned int lingno2_membership; 

float smallest; 
float smallest2; 

rule_node - 0; 

for (lingnol_membership - 0; lingnol_membership < MEMBERSHIP_N01; 
lingno 1 _membership++) 

{ 
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smallest *» ActJ2[inputJingui_nol]|lingnol .membership]; 


for (lingno2_membership - 0; lingno2_membership < MEMBERSHEP.N02; 
lingno2_membership++) 

{ 


/* printf(“rule_node -%d\t”,rule_node); */ 

smallest2 - Act_12[input_lingui_no2][lingno2_membership]; 

/* printf(“smallest smallest); */ 

/* printf(“smallest2 smallest2); */ 

if (smallest < smallest2) 

{ 

Act_13[rule_node] - smallest; 

/* printf(“Act_13 **%fSn”, Act_13[rule_node]); */ 

> 

else 

{ 

Act_13[rule_node] - smallest2; 

/* printf(“Act_13 -%f\n”, Act_13[rule_node]); */ 

> 

rule_node++; 

if (rule.node > RULE_NODE) 

{ 

printf(“ERROR number of rule nodes are greater than the expected value\n”); 
return (FALSE); 

> 

/* getchQ; */ 


> /* End of lingno2_membership */ 

) /* End of lingno 1 .membership */ 

return (TRUE); 

} /* End of the program */ 
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#include <stdio.h> 

#include <st<^lib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

/******************************************************************************* 

* function get_act40 

* 

* This function calculates the output of the layer two of the network. This 

* function uses a function pointers points to function activation_12. 

* 

* Input parameter none. 

* 

* Return value: int type. 

* 

*programmerLeung Kam LungMS. Thesis9/28/93 
* 

*/ 


int 

get_act4(unsigned int output_lingui_no, unsigned int index) 

{ 

int retval; 

int (*act_function_ptr) (float, unsigned int); 
act_function_ptr - activation_14; 

retval - (*act_function_ptr) (Output[outputJinguijno][index], output_lingui_no); 
if (retval — FALSE) 

< 

printf(“reture value of act_function_ptr is > or < max,min doubleNn”); 
return (FALSE); 

> /* End of if retval */ 

return (TRUE); 

> 

/ 

*************. & a;******************************************************************* 
*************** 

* This function is called by get_act2 to set up the parameters for calculates the output value 


136 



* for layer of the network. 

* 

' A 

* Input parameter inputx, 

*input_lingui_no. 

* 

* Return value: TRUE if successful! calculate the output values. 

* FALSE otherwise. 

* 

* 

♦programmer.Leung Kam LungMS. Thesis9/28/93 
*/ 

int 

activation_14(float output, unsigned int output_lingui_no) 

{ 

int membership; 
float temp4; 

double (*activation_14_ptr) (float, float, float); 
activation_14_ptr = 14_membership; 

for (membership = 0; membership < Output_membership_array[output_lingui_no]; membership++) 

{ 


temp4 - Act_14[output_lingui_no] [membership] - (float) (*activation_14_ptr) 
(Mean_output[output_lingui_no][membership], 
Output_term_Deviation[output_lingui_no][membership], 
output); 

/* printf(“act4 %f\n”,temp4); *1 
if (temp4 > MAXFLOAT) 

{ 


printf(“Act_14 value of output_lingui_no %d membership %d is > or < expected value\n”, 
output_lingui_no, membership); 

return (FALSE); 

> 


> /* End of for membership */ 

return (TRUE); 
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> 


/* End of activation_12 */ 


/****************************************************************************** 


* The function calculates the activation value for membership function in 

* layer 4 of the NNFLC. 

* input parameters: 

♦mean 
♦deviation 
*u(input value) 

* output: 

♦Return type double 


* 


*programmer:Leung Kam LungMS. Thesis9/28/93 
* 


*/ 


double 

14_membership(float mean, float deviation, float u) 

{ 


/* 

* printf(“4 deviation %f, mean %f, float %f”, deviation, mean, u); 

* getchar(); 

*/ 

return (exp(-pow(((double) u - (double) mean), 2) / (pow((double) deviation, 2)))); 


> 
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#include <stdio,h> 

#include <stcflib.h> 

#include <stdio.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file'^Iobalex.h” 


int 

find_max_w4(void) 

{ 

float weight4, 
max_w4; 

unsigned int output_lingui_no, 
output_term_membership, 
rule_index, 
max_w4_index; 


for (output_lingui_no *~ 0; output_lingui_no < OUTPUT_LINGUI_NO; outputJingui_no++) 
{ 


for (rule_index - 0; rule_index < RULE_NODE; rule_index++) 

{ 


/* 

* printf(“output_term_membership %d, rule 

* %d”,output_term_merabership, rule_index); 

*/ 

/* 

* printf(“ weight4 %f 

* “,W J4[output_lingui_no] [output_term_membership][rule_index]); 

*/ 


max_w4 - -99999.0; 

for (output_term_membership - 0; output_term_membership < 
Output_membership_anay[output_lingui_no]; 
output_term_membership++) 

{ 


weight4 - WJ4[output_lingui_no][output_term_membership][rule_index]; 



/* printf(“weight4 %f\n”,weight4); */ 

*'J% 

if (weight4 >- max_w4) 

{ 

max_w4 - weight4; 

max_w4_index - output_term_membership; 

/* 

* printf(“ max_w4_index %d, w_4 %f\n “,max_w4_index, 

* weight4); 

*/ 


if (weight4 > MAXFLOAT) 

{ 

printf(“weight4 value of output_lingui_no %d output term membership %d rule_index &d 
is > or < expected value\n”, 

output_lingui_no, output_term_membership, rule_index); 
return (FALSE); 


> 


> /* End for output_term_membership */ 


for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 
output_term_membership++) 

{ 

if (output_term_membership — - max_w4_index) 

{ 

W_14[output_lingui_no][output_term_membership][rule_index] - 1.0; 

> 

else 

{ 

WJ4[output_lingui_no][output_term_membership][rule_index] - 0.0; 

} 
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> 


} . , /* End second for output_term_membership */ 

/* End rule_index */ 

} /* End for output_Iingui_no */ 

return (TRUE); 

} /* End of function change_w4 */ 
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#include <stdio.h> 

#include <st(ilib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalexJi” 

^*******^** 5 f: 3 f:* + + **** + + + + **** ;1 c + + + *******^*****************:***** + *** ************* 

* function find_output() 

* 

* This function calculates the output of the layer two of the network. This 

* function uses a function pointers points to function activation_12. 

* 

* Input parameter none. 

* 

* Return value: int type. 

* 

♦programmer.Leung Kam LungMS. Thesis9/28/93 
* 

*/ 


int 

findoutput(void) 

{ 

int retval2, 
retval3, 
retval4; 

unsigned int input_lingui_no, 
input_lingui_nol , 
input_lingui_no2, 
output_lingui_no; 

float outl, 
f, 
f2; 

unsigned int index, 
output_term_membership, 
membership; 

int (*act_function_ptr2) (unsigned int, unsigned int); 

int (*act_function_ptr3) (unsigned int, unsigned int); 

int ( * act_function_downup_ptr4) (void); 
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input_lingui__nol = 0; 
input_Iingui_no2 - 1 ; 
output_lingui.no - 0; 


act_function_downup_ptr4 - get_downup_act4; 
act_function_ptr3 = get_act3; 
act_function_ptr2 - get_act2; 

for (index - 0; index < MAX_INPUT_INDEX; index++) 

{ 


for (input_lingui_no = 0; input_lingui_no < INPUT_LINGUI_NO; input_lingui_no++) - 

< 


retval2 » (*act_function_ptr2) (input_lingui_no, index); 

for (membership - 0; membership < Input_membership_array[input_lingui_no]; member- 
ship++) 

{ 

fprintf(Act2_fptr, “%1.6f Act_12[input_lingui_no] [membership]); 

> 


> /* End of for input_lingui_no */ 
fprintf(Act2_fptr, “\n”); 

retval3 - (*act_function_ptr3) (input_lingui_no 1 , input_lingui_no2); 
for (membership - 0; membership < RULE_N0DE; membership++) 

< 

fprintf(Act3_fptr, “%1.6f Act_13 [membership]); 

> 

fprintf(Act3_fptr, ‘V’’); 

retval4 - (*act_fimction_downup_ptr4) (); 

for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


for (membership - 0; membership < Output_membership_array[output_lingui_no]; member- 
ship++) 

{ 

fprintf(Act4_fptr, “%1.6f Act downup_14[output_lingui_no][membership]); 

> 


> 

fprintf(Act4_fptr, “\n”); 
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for (output_lmgui_.no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


output_down_up[index][output_lingui_no] - 0.0; 
f - f2 - 0.0; 

for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 
output_term_membership++) 

{ 


f +- (Mean_output[output_lingui_no][output_tenn_membership] * 

Output_term_Deviation[output_lingui_no] [output_term_membership] * 
Act_downup_14[output_lingui_no][output_term_membership]); 

f2 +- (Output_term_Deviation[output_lingui_no][output_term_membership] * 
Act_downup_14[output_lingui_no][output_term_membership]); 

> /* End of output_term_membership */ 

outl - output_down_up[index][output_lingui_no] - f / f2; 

/* 

* printf(“ actual output %f Output[0][index]); printf(“ the 

* calculated output is %f\n”, outl); getchar(); 

*/ 

> /* End of output_lingui_no */ 

if (retval2 — FALSE) 

{ 

printf(“reture value of act_functjon_ptr2 is > or < max,min doubleNn”); 
return (FALSE); 

> /* End of if retval */ 

if (retval3 — FALSE) 

{ 

printf(“reture value of act_function_ptr3 is > or < max,min double\n”); 
return (FALSE); 

} /* End of if retval */ 

if (retval4 — FALSE) 

{ 

printf(“reture value of act_function_downup_ptr4 is > or < max,min double\n”); 
return (FALSE); 

} /* End of if retval */ 
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> 


/* End of for input_index */ 


* i» 

fclose(Act2_fptr); 

fclose(Act3_fptr); 

fclose(Act4_fptr); 

return (TRUE); 

} 
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#include <stdio.h> 

#include <stdlib.h> 

#include “c:\borlandc\file\thesis.h” 
#inchide“c:\borlandc\file\globalex.h” 









j ; 

I 

U 


I! 


I***************************************************************************** 


* function init_phasel() 

* 


* This function excute all initial procedure for getting the center of the 

* membership and the deviation of input and output term nodes.So after this 

* function the input signal can be pass through the membership on layer3 and 

* layer4. 

* 


♦Input parameter inputxl, 

*inputx2. 

* 

* Return value: TRUE if each of the center of membership and the deviation 
♦are calculated. 

♦FALSE otherwise. 

* 

* Programmer 

* Leung Kam Lung 9/30/93MS. Thesis. 

* 


*/ 


int 

initjphasel(void) 


inti, 

ind; 

float time, 
time2; 

unsigned int input_lingui_no, 
output_lingui_no; 

ind - 1; 
time - 0.10; 
time2 - 0.10; 
initials_all(); 
mput_lingui_.no - 0; 
output_lingui_.no - 0; 

while ((time > 0.000001) II (time2 > 0.000001)) 

{ 

for (i - 0; i < MAX_INPUT_INDEX; i++) 
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{ 

updaje_input_12_w(Inputx[input_lingui_no][i], time, input_lingui_no); 
update_input_12_w(Inputx[input_lingui_no + l][i], time, input_lingui_no + 1); 
update_output_w(Output[outputJingui_no][i], time2, output_Iingui_no); 

} 

time ” time * 0.95; 
time2 - time2 * 0.95; 


ind++; 

> 

printf(“ ind %d\n”, ind); 

if (!save_mean_input()) 

{ 

printf(“Mean_input FALSEXn”); 

> 

if (!save_mean_output(» 

{ 

printf(“Mean_output FALSENn”); 

> 


I* 

* up to now only the center(mean) of the input and output membership are 

* calculated 
*/ 

for (input_lingui_no - 0; input_lingui_no < INPUT JJNGUI_NO; input_lingui_no++) 

{ 


if (find_input_deviation(input_lingui_no)) 

{ 

printf(“ Deviation of each Membership is find for INPUT_TERM_NODE %dsn”, 
input_lingui_no); 

> 

else 

{ 

printf(“ NOT Deviation of each Membership is find for INPUT_TERM_NODE %dSn”, 
input_lingui_no); 

return (FALSE); 

> 
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} . /* End for input_lingio_no */ 

if (!save_input_deviation()) 

{ 

printf(“CAN NOT SAVE INPUT DEVIATION^”); 

> 

for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


if (find_output_deviation(output_lingui_no)) 

< 

printf(“ Deviation of each Membership is find for OUTPUT_TERM_NODE %d\n”, 
output_lingui_no); 

> 

else 

{ 

printf(“ NOT Deviation of each Membership is find for OUTPUT_TERM_NODE %d\n”, 
output_lingui_no); 

return (FALSE); 

> 

> /* End for output_lingio_no */ 

if (!save_output_deviation()) 

{ 

printf(“CAN NOT SAVE OUTPUT DEVIATION^”); 

> 

/* 

* now all deviation of each membership for input term nodes and output 

* term nodes are calculated 
*/ 


if (!update_w_14()) 

< 

printf(“udate_w_14 FALSEsn”); 
return (FALSE); 

} 

if (!save_w4()) 

{ 

printf(“can not save weight 4\n”); 
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return (FALSE); 

> 

if (!find_output()) 

{ 

printf(“ Can not find the output'll”); 
return (FALSE); 

> 


return (TRUE); 



#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include “c:\borlandc\file\thesis.h” 
#include“c:\borlandc\file\globalex.h” 


int 

save_w4(void) 

{ 

float weight4; 

unsigned int output_lingui_no; 
unsigned int output_term_membership; 
unsigned int rule_index; 

for (output _lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 

for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 
output_teim_membership++) 

{ 


for (rule_index = 0; rule_index < RULE_NODE; rule_index++) 

{ 

fprintf(W_14_fptr, “%1.6f 

W_14[output_lingui_no] [output_term_membership] [rule_index]); 

} /* End rule_index */ 

fprintf(W_14_fptr, “\n”); 

> /* End for output_term_membership */ 

> I* End for output_lingui_no *1 

fclose(W_14_fptr); 

return (TRUE); 

) /* End of function change_w4 */ 
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#include <stdio.h> 

#include <stdlib.h> 

#include <niath.h> 

#include “c:\borlandc\file\thesis.h” 

#include “c:\borlandc\file\globalex.h” 

int 

save_w4(void) 

{ 

float weight4; 

unsigned int output_lingui_no; 
unsigned int output_term_membership; 
unsigned int rule_index; 

for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


for (output_term_membership - 0; output_term_membership < 
Output_membership_array[output_lingui_no]; 
output_term_membership++) 

{ 

. for (rule_index - 0; rule_index < RULE_NODE; rule_index++) 

{ 


fprintf(W_14_fptr, "%d (int) 

W_14[output_lingui_no][output_term_membership][rule_index]); 

> /* End rule_index */ 

fprintf(W_14_fptr, “\n”); 

} /* End for output_term_membership */ 

> /* End for output_lingui_.no */ 

if (fclose(W_14_fptr)) 

{ 

return (FALSE); 

> 

return (TRUE); 


} /* End of function change_w4 */ 
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#include <stdio.h> 

#include <st$lib.h> 

#include “c:\borlandcM51e\ihesis.h” 
#include “c:\borlandc\file\globalex.h” 


int 

update_w(int prt_index) 

{ 

unsigned int centerjndex, 
devi_index, 
input_lingui_no, 
output_lingui_no, 
membership; 
int go_print; 

if (prt_index = 249) 

{ 

go_print - TRUE; 

> 

else 

{ 

go_print =■ FALSE; 

> 


for (output_lingui_no - 0; output_lingui_no < OUTPUT_LINGUI_NO; output_lingui_no++) 

{ 


for (center_index - 0; center_index < Output_membership_array[output_lingui_no]; 
center_index++) 

{ 

Mean_output[output_lingui_no] [center_index] - Mean_output[output_lingui_no] [center_index] 

Mean_output_delt[output_lingui_no][center_index]; 

if (go_print — TRUE) 

{ 

fprintf(Final_mean_output_fptr, “%2.6f Mean_output[outputJingui_no] [centerjndex]); 

> 

Mean_output_delt[output_lingui_no][center_index] - 0.0; 

} /* center_index */ 


152 



if (go_grint *~= TRUE) 

{ 

fprintf(Final_mean_output_fjptr, “\n”); 

> 

for (devi_index - 0; devi_index < Output_membership_array[output_lingui_no]; 
devi_index++) 

{ 


Output_term_Deviation[output_lingui_no] [devi_index] - 
Output_term_Deviation[output_lingui_no] [devijndex] + 

Output_term_Deviation_curr[output_lingui_no]tdevi_index]; 

if (go_print — « TRUE) 

{ 

fprintf(Final_output_tenn_deviation_fptr, “%1 .6f 
Output J:ermJDeviation[outputJingui_no][devi_index]); 

> 

Output_term_Deviation_curr[output_lingui_no][devi_index] - 0.0; 

) /* End for devi_index */ 

if (go_print — TRUE) 

{ 

fprintf(Final_output_term_deviation_fptr, “\n”); 

> 


> /* End for output_lingui_no */ 

for (input_lingui_no -= 0; input_lingui_no < INPUT_LINGUI_NO; input_lingui_no++) 

< 


for (membership - 0; membership < Input_membership_array[input_lingui_no]; 
membership++) 

{ 


/* 

* printf(“Mean_input no %d mem_ship %d %2.5f 

* “,input_lingui_no,membership, 

* Mean_input[input_lingui_no] [membership]); 

*/ 

Mean_input[input_lingui_no)[membership] *» Mean_input[input_lingui_no][membership] + 
Mean_input_curr[input_lingui_no][membership]; 
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if (go_print — - TRUE) 

< 

fprintf(Final_mean_input_fptr, “%2.6f Mean_input[input_lingui_no] [membership]); 

> 

/* 

* printf(“delta %2.5f 

* Mean_input_curr[input_lingui_iio][membership]); 

*1 

/* 

* printf(“FMean_input %2.5f \n”, 

* Mean_input[input_lingui_no][membership]); 

*/ 

Mean_input_curr[input_lingui_no] [membership] - 0.0; 

Input_term_Deviation[input_lingui_no] [membership] - 
Input_term_Deviation[input_lingui_no] [membership] + 

Input_term_Deviation_curr[input_lingui_no][membership]; 


if (go_print — TRUE) 

{ 

fprintf(Final_input_term_deviation_fptr, “% 1 .6f 
Input_term_Deviation[input_lingui_no][membership]); 

> 

/* 

* printf(“delta Mean_Devi input_l_no %d mem_ship %d 

* %2.5f\n”,input_lingui_no,membership, 

* Input_term_Deviation_curr[input_lingui_no][membership]); 

*/ 

Input_term_Deviation_curr[input_lingui_no][membership] - 0.0; 

> /* End for membership */ 

} /* End for input_lingui_no */ 

if (go_print — TRUE) 

{ 

fprintf(Final_mean_input_fptr, “\n”); 
fprintf(Final_input_tenn_deviation fptr, * V’); 

> 

return (TRUE); 
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> 


/* End of update_w */ 
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#include <stdio.h> 

#include <st<Jlib.h> 

#include <time.h> 

#include “c:'^)orlandc\fUe\thesis .h” 
#include “c:\borlandc\file\globalex.h” 


main(void) 


< 


unsigned int index, 
input_index, 
find, 
counter, 
total; 

float input 1, 
input2, 
outputl, 
max_inputl, 
max_input2, 
max_outputl, 
min_inputl, 
min_input2, 
min_outputl, 
input_slopel, 
input_slope2, 
output_slopel; 

float number, 

long big, 
small; 

float array_input 1 [M AX_INPUT_INDEX] , 
array _input2 [M AX_INPUT_INDEX] , 
array_outputl [MAX_INPUT_INDEX] ; 

float array _counter[MAX_INPUT_INDEX]; 

FILE *Error_in; 

FILE *Fback_in; 

FILE *Consigna_out; 

FILE *Max_input_value_fptr, 

FILE *Slope_fptr; 
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if ((Error Jn - fopen(“error.dat”, “r”)) — NULL) 

< 

fprintf(stderr, “CAN NOT OPEN ERROR FILEmi”); 
return (FALSE); 

> 

if ((Fbackjn - fopen(“fback.dat”, “r”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN FBACK FILEMi”); 
return (FALSE); 

> 

if ((Consigna_out - fopen(“consigna.dat”, “r”» —> NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN CONSIGNA FILEVn”); 
return (FALSE); 

> 


for (index - 0; index < MAX JNPUT JNDEX; index++) 

{ 

array_counter!index] - -10.0; 

> 


total - 0; 


for (index - 0; index < MAXJNPUTJNDEX; index++) 

{ 


fscanf(Error_in, “%f’, &inputl); 
fscanf(Fback_in, “%f &input2); 
fscanf(Consigna_out, “%f\ &outputl); 

do 

{ 

find -FALSE; 

counter = random(MAXJNPUTJNDEX); 
if (array_counter[counter] -- -10.0) 

{ 

array_counter[counter] - 10.0; 
array_inputl [counter] - inputl; 
array_input2[counter] - input2; 
airay_outputl [counter] -outputl; 
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find = TRUE; 

* .1 

total++; 

> 

> while (find — FALSE); 


> 


/********************************************************** ************ 


* This next two for loop is to normalize the trainning data to a range 
*of-l and 1. 


* 


*/ 


maxjnputl - array_inputl[0]; 
max_input2 - array_input2[0j;; 
max_outputl - array_outputl [0]; 

min_inputl - arrayjnputl[0]; 
min_input2 - array_input2[0]; 
min_outputl - array_outputl [0]; 

for (index - 1; index < MAX_INPUT_INDEX; index++) 

{ 


if (maxjnputl < array Jnputl [index]) 

{ 

maxjnputl - array Jnputl [index]; 

} 

if (minjnputl > array Jnputl [index]) 

{ 

minjnputl - array Jnputl [index]; 

> 


if (maxjnput2 < array Jnput2[index]) 

{ 

maxjnput2 - array Jnput2[index]; 

> 

if (minjnput2 > array Jnput2[index]) 

{ 

minjnput2 - array Jnput2[index]; 

> 
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if (max outputl < array _outputl [index]) 

{ 

max_outputl « array_outputl [index]; 

> 

if (min_outputl > array_outputl [index]) 

{ 

min_outputl - array_outputl [index]; 

> 


> /* End for index loop */ 

printf(“max_outputl %f max_inputl %f max_input2 %f\n”, max_outputl, max_inputl, 
max_input2); 

printf(“min_output 1 %f min_inputl %f min_input2 %f\n”, min_outputl, min_inputl, min_input2) 
getchar(); 

if (fabs(max_inputl) > fabs(min_inputl)) 

{ 

input_slopel - 1.0/fabs(max_inputl); 

> 

else 

{ 

input_slopel - 1.0/ fabs(min_inputl); 

> 

if (fabs(max_input2) > fabs(min_input2)) 

{ 

input_slope2 -1.0/ fabs(max_input2); 

> 

else 

{ 

input_slope2 - 1.0/ fabs(min_input2); 

> 

if (fabs(max_outputl) > fabs(min_outputl)) 

{ 

output_slopel - 1.0 / (fabs(max_outputl *2.0)); 

} 

else 

< 

output_slopel » 1.0 / (int) ((fabs(min_outputl * 2.0)) + 1); 

> 
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for (index,” 0; index < MAX JNPUTJNDEX; index++) 

{ 


arrayjnputl [index] - array_inputl [index] * input_slopel; 
array_input2(index] - array_input2[index] * input_slope2; 
array_outputl [index] - 0.5 + (array_outputl [index] * output_slopel); 


> 

max_inputl - array_inputl[0]; 
max_input2 - axray_input2[0];; 
max_outputl » array _outputl[0]; 

min_inputl - array_inputl[0]; 
min_input2 - array_input2[0]; 
min_outputl » array_outputl[0]; 

for (index = 1; index < MAX_INPUT_INDEX; index++) 

{ 


if (max_inputl < array_inputl [index]) 

{ 

max_inputl - arrayjnputl [index]; 

> 

if (min_inputl > airay_inputl [index]) 

{ 

min_inputl - arrayjnputl [index]; 

> 


if (max_input2 < array_input2[index]) 

{ 

max_input2 - array _input2[index]; 

} 

if (min_input2 > array_input2[index]) 

{ 

min_input2 - array_input2[index]; 

> 


if (max_outputl < array_outputl [index]) 

{ 

max_outputl - array_outputl [index]; 
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> 


if (mirijDutputl > array_outputl [index]) 

{ 

min_outputl - array_outputl [index]; 

> 


} /* End for index loop */ 

printf(“max_outputl %f max_inputl %f max_input2 %f\n”, max_outputl, raax_inputl, 
max_input2); 

printf(“min_outputl %f min_inputl %f min_input2 %f\n”, min_outputl, min_inputl, min_input2); 
getchar(); 

fclose(Error_in); 

fclose(Fback_in); 

fclose(Consigna_out); 


if ((Max_input_value_fptr - fopen(“maxinput.txt”, “w”)) == NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN MAX INPUT VALUE FELE\n”); 
return (FALSE); 

> 


if ((Error_fptr -* fopen(“error.txt”, “w”)) — ■ NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN ERROR FILE\n”); 
return (FALSE); 

> 

if ((Fbackjfptr - fopen(“fback.txt”, “w”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN FBACK FILE\n”); 
return (FALSE); 

> 

if ((Consigna_fptr - fopen(“consigna.txt”, “w”)) — NULL) 

< 

fprintf(stderr, “CAN NOT OPEN CONSIGNA FILE\n”); 
return (FALSE); 

} 

if ((Slopejptr - fopen(“slope.txt”, “w”)) — NULL) 
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{ 

fprintf(stderr, “CAN NOT OPEN CONSIGNA HLEmT); 
return (IrALSE); 

> 


fprintf(Max_input_value_fptr, “% 2.7e\n”, max_inputl); 
fprintf(Max_input_value_f^tr, “% 2.7e\n”, min_inputl); 
fprintf(Max_input_value_fptr, “% 2.7 eW\ max_input2); 
fprintf(Max_input_value_fptr, “% 2.7e\n”, min_input2); 
fprintf(Max_input_value_fptr, “% 2.7eW, max_outputl); 
fprintf(Max_input_value_f^tr, “% 2.7e\n”, min_outputl); 

for (index - 0; index < MAX_INPUT_INDEX; index++) 

< 


number - array_inputl [index]; 
small - ((long) (number * 100000) * 10); 
big - (long) (number * 1000000); 
if ((big - small) >= 5) 

{ 

number - (small +- 10); 

> 

else 

{ 

number - small; 

> 

/* fprintf(Error_fptr,”% 2.7e\n”,array_inputl [index]); */ 
fprintf(Error_fptr, “% 2.7f\n”, number / 1000000.0); 

number - array_input2[index]; 
small - (Gong) (number * 100000) * 10); 
big - Gong) (number * 1000000); 
if ((big - small) >- 5) 

{ 

number - (small +- 10); 

> 

else 

{ 

number - small; 

} 

/* fprintf(Fback_fptr,”% 2.7eV’,array_input2[index]); */ 
fprintf(Fback_fptr, “% 2.7f\n”, number / 1000000.0); 
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number - array _outputl [index]; 


small “’((long) (number * 100000) * 10); 
big “ (long) (number * 1000000); 

if ((big - small) >- 5) 

< 

number - (small +« 10); 

> 

else 

{ 

number - small; 

> 

/* fprintf(Consigna_fptr,”% 2.7e\n”,array_outputl [index]); */ 
fprintf(Consigna_fptr, “% 2.7f\n”, number / 1 000000.0); 


> 


fprintf(Slope_fptr, “% 2.7f\n”, input_slopel); 
fprintf(Slope_fptr, “% 2.7fm”, input_slope2); 
fprintf(S lope_fptr, “% 2.7f\n”, output_slopel); 

fclose(Error_fptr); 

fclose(Fback_fptr); 

fclose(Consigna_fptr); 

fclose(Max_input_value_fptr); 

fclose(Slope_fptr); 

return (TRUE); 

> 
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#include <stdio.h> 

#include “c:\borlandc\file\thesis.h” 
#include “c:\borlandc\file\globalex-h” 


int 

read_input(void) 

{ 

int inputjndex; 

for (input_index - 0; inputjndex < MAX JNPUTJNDEX; input Jndex++) 

{ 


fscanf(Error_fptr, “%f’, &lnputx[0][inputjndex]); 
fscanf(Fback_fptr, “%f\ &Inputx[l][inputJndex]); 
fscanf(Consigna_fptr, “%f’, &Output[0] [inputjndex]); 

/* 

* printf(“Error_ptr %f\n”,lnputx[0][inputjndex]); printf(“Fback_ptr 

* %f\n”,Inputx[l][input_index]); printf(“Consigna_ptr 

* %f’,Output[0][input_index]); getchar(); 

*/ 


fclose(Error_fptr); 

fclose(Fback_fptr); 

fclose(Consigna_fptr); 

if ((Slope Jptr = fopen(“slope.txt”, “r”)) — NULL) 

{ 

fprintf(stderr, “CAN NOT OPEN SLOPE FILENn”); 
return (FALSE); 

> 


fscanf(Slope_fptr, “%f\ &Input_slopel); 
fscanf(SlopeJptr, “%f &Input_slope2); 
fscanf(Slope_fptr, “%f\ &Output_slopel); 
fclose(Slope_fptr); 

if ((Max_data_fptr - fopen(“maxinputtxt”, “r”)) — NULL) 
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{ 

fprintf(stderr, “CAN NOT OPEN MAX INPUT FBLE\n”); 
return (FALSE); 

} 


for (inputjndex - 0; inputjndex < ((2 * INPUT_LINGUI_NO) + (2 * OUTPUT_LINGUI_NO)); 
input_index++) 

{ 

fscanf(Max_data_fptr, “%f &Max_inputl); 
fscanf(Max_data_fptr, “%f’, &Min_inputl); 
fscanf(Max_data_fptr, “%f &Max_input2); 
fscanf(Max_data_f^tr, “%f &Min_input2); 
fscanf(Max_data_fptr, “%f\ &Max_oulputl); 
fscanf(Max data_fptr, “%f\ &Min_outputI); 

> 

printf(“Max_inputl %f \n”, Max_inputl); 
printf(“Min_inputl %f\n”, Min_inputl); 
printf(“Max_input2 %f Nn”, Max_input2); 
printf(“Min_input2 %f \n”, Min_input2); 
printf(“Max_outputl %f\n”, Max_outputl); 
printf(“Min_outputl %f \n”, Min_outputl); 

fclose(Max_data_fptr); 

getcharQ; 

return (TRUE); 

> 
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